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Overview 


This  final  report  summarizes  the  theoretical  and  experimental  activities 
carried  out  in  joint  collaboration  by  the  Quantum  Optics  Group  at  Drexel 
University  and  the  Nonlinear  Dynamics  and  Laser  Physics  Group  at  Bryn  Mawr 
College  during  the  period  October  15,  1985  -  August  31,  1987.  Although  much  of 
the  research  conducted  under  this  constract  has  involved  a  close  interaction  between 
theorists  and  experimentalists,  for  convenience  we  list  the  main  accomplishments 
under  separate  headings  of  theory,  experiments  and  special  activities. 

The  following  summaries  represent  highlights  of  technical  results  which  are 
discussed  in  greater  depth  in  the  enclosed  publications  by  the  members  of  our 
groups. 


The  Role  of  Transverse  Effects  in  Laser  Instabilities 

1.  Introduction 

One  of  the  unexpected  features  of  early  laser  systems  was  the  appearance  of  output 
pulsations  even  under  steady  or  nearly  steady  pumping  conditions.  In  fact,  spiking  had  already 
been  observed  in  maser  experiments*  even  before  the  discovery  of  the  laser,  with  the  proliferation 
of  solid  state  optical  devices,  this  effect  acquired  the  status  of  a  nearly  universal  feature^.  The 
interpretation  of  these  phenomena  as  direct  manifestations  of  intrinsic  dynamic  instabilities,  linked 
to  the  nonlinear  character  of  the  emission  should  be  ranked  among  the  most  significant  theoretical 
advances  of  the  late  1950's  and  early  1960's^’^. 

The  most  common  modem  framework  for  the  description  of  laser  dynamics  is  the 
plane-wave  Maxwell-BIoch  model,  a  formulation  flexible  enough  to  conform  to  most  experimental 
conditions  and  capable  of  yielding  a  detailed  insight  into  the  laser  working  mechanism.  Yet,  for  all 
the  important  contributions  that  this  model  has  made  to  laser  physics,  it  has  also  shown  a  consistent 
pattern  of  quantitative  and  even  qualitative  disagreement  with  experimental  facts,  especially  when 
the  instabilities  become  a  dominant  dynamical  feature^. 

An  important  starting  point  for  our  discussion  is  the  recognition  that  the  plane-wave 
Maxwell-BIoch  equations  do  predict  the  existence  of  instabilities  leading  to  both  regular  and  chaotic 
pulsations;  this  is  a  well  established  theoretical  fact^.  Thus,  homogeneously  broadened  systems 
are  known  to  produce  output  oscillations  when  operating  both  in  single  and  multimode 
configurations,  and  in  and  out  of  resonance,  relative  to  the  center  of  the  atomic  gain  line.  Of  special 
significance  is  the  existence  of  an  isomorphism  between  the  single-mode  laser  model  and  the 
Lorenz  equations^3.  The  Lorenz  equations,  originally  derived  to  simulate  the  onset  of  convective 
instabilities^,  are  well  known  in  the  mathematical  literature  as  the  paradigm  of  deterministic  chaos^. 
Thus,  the  isomorphism  between  the  laser  and  the  Lorenz  equations  supports  the  notion  that  the 
laser  itself  can  be  the  source  of  chaotic  behavior. 


While  the  link  between  the  observed  laser  pulsations  and  dynamic  instabilities  is  a  very 
reasonable  theoretical  proposition,  the  prediction  that  these  phenomena  can  occur  only  for  very 
large  values  of  the  pump  parameter,  especially  under  resonant  conditions,  is  the  source  of  major 
difficulties  with  the  interpretation  of  experimental  results.  Indeed,  if  one  took  the  theoretical 
threshold  values  too  seriously,  one  would  have  to  conclude  that  laser  instabilities  are  observable 
only  under  exceptional  circumstances. 

An  additional  difficulty  with  the  plane-wave  Maxwell-Bloch  equations  is  the  unconditional 
stability  of  their  rate  equations  limit.  Accepting  this  result  leads  to  the  conclusion  that  all  laser 
instabilities  are  manifestations  of  atomic  coherence,  this  assertion  is  hardly  in  agreement  with  the 
behavior  of  Ruby,  Nd:YAG,  C02  and  other  lasers  for  which  the  validity  of  the  rate  equations 
description  is  well  justified. 

Better  agreement  between  theory  and  experiments  is  obtained  when  the  laser  undergoes  a 
phase  instability^,  a  phenomenon  that  requires  operation  under  detuned  conditions  and  which  is 
responsible,  for  example,  for  the  observed  mode  hopping  in  C02  laser  systems^.  However,  the 
physical  mechanism  that  drives  a  phase  instability  is  quite  different  from  that  which  is  operative 
under  resonant  conditions  (i.e.  the  amplitude  instability) 111  so  that  the  theory  of  phase  instabilities  is 
not  very  helpful  in  clarifying  the  nature  of  the  "high  second  threshold  problem".  The  theory  of 
inhomogeneously  broadened  lasers  is  also  in  good  qualitative  agreement  with  the  observed  unstable 
phenomena1 in  this  case,  the  lowering  of  the  second  laser  threshold  to  the  observed  experimental 
range1  ^  is  probably  due  to  the  large  increase  in  the  number  of  degrees  of  freedom  that  are  needed  to 
describe  the  laser  dynamics,  and  offers  little  clue  as  to  the  failing  of  the  homogeneously  broadened 
model. 

It  is  probably  fair  to  say  that  the  plane- wave  approximation,  long  viewed  as  adequate  in 
capturing  the  essential  aspects  of  laser  dynamics,  may  be  a  much  stronger  theoretical  constraint 
than  anticipated.  In  fact,  important  warning  signals  in  this  direction  came  from  the  lack  of 
quantitative  agreement  between  the  predictions  of  the  plane-wave  studies  of  optical  bistability1  ^  and 
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the  careful  absolute  measurements  performed  by  Kimble  and  collaborators  under  steady  state 
conditions^.  If,  however,  these  theoretical  propositions  were  in  disagreement  with  the 
experiments  in  a  quantitative  way,  more  serious  qualitative  discrepancies  have  emerged  from  the 
study  of  self-pulsing  in  bistable  systems.  Here  the  plane-wave  theory  predicts  a  rich  variety  of 
dynamical  behaviors ^  which  do  not  appear  to  have  a  counterpart  among  the  observed  patterns^. 

The  main  objective  of  this  research  is  to  produce  evidence  that  transverse  degrees  of  freedom 
play  a  far  more  important  role  in  laser  dynamics  than  might  be  expected.  In  fact,  we  have  known 
for  some  time  that  transverse  effects  can  have  a  strong  influence  on  the  stationary  and  dynamic 
behavior  of  passive  driven  systems^  and  that  plane-wave  and  transverse  models  can  approach 
chaos  by  way  of  very  different  routes^.  The  combined  message  of  a  small  but  growing  number 
of  experimental^  and  theoretical  reports^  supports  the  belief  that  effects  related  to  a  departure 
from  the  plane-wave  configuration  may  be  even  more  influential  when  the  optical  resonator 
contains  an  active  rather  than  a  passive  medium.  If  we  accept  this  premise,  the  inclusion  of 
transverse  degrees  of  freedom  into  the  existing  theories  becomes  a  necessary  requirement  for  a 
successful  description  of  realistic  laser  systems. 

Indeed,  several  theoretical  studies  have  already  attempted  to  improve  on  the  plane-wave 
description.  An  analysis  of  the  modes  and  of  their  spatial  stability  in  a  realistic  resonator  can  be 
found,  for  example,  in  Ref.  21.  Reference  22  contains  a  study  of  the  emission  frequency  of  an 
inhomogeneously  broadened  laser  when  the  field  configuration  is  no  longer  of  the  pure  Gaussian 
type.  The  problem  of  laser  dynamics  in  the  presence  of  transverse  effects  was  addressed  in  Ref.  23 
using  an  approximate  set  of  equations  and  evidence  was  presented  for  the  appearance  of 
self-pulsing. 

A  useful  alternative  approach  is  based  on  the  assumption  that  the  cavity  field  maintains  a 
fixed  Gaussian  profile  during  the  evolution  17b, 24  On  the  surface  this  appears  to  be  a  rather  crude 
approximation  because  the  presence  of  an  active  medium  may  be  expected  to  introduce  significant 
distortion  in  the  transverse  pattern  of  the  field.  In  the  case  of  passive  systems,  however,  this 
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approach  has  yielded  quite  reasonable  results^,  perhaps  because  of  the  constraints  that  the  injected 
mode-matched  signal  imposes  on  the  cavity  field.  The  situation  is  likely  to  be  more  critical  with 
active  systems,  so  that  a  more  realistic  approach  may  be  required.  In  fact,  even  the  simple 
Gaussian  model  of  the  laser  yields  remarkable  deviations  from  the  plane- wave  behavior,  the  Risken 
and  Nummedal  instability^’0  of  the  longitudinal  sidebands,  for  example,  is  suppressed  under 
resonant  conditions  if  the  cavity  field  has  a  TEMqq  profiled 

In  this  study  we  consider  those  physical  and  geometrical  features  that  are  absent  by  definition 
from  the  plane-wave  theory.  To  be  specific,  we  focus  on  the  effects  of  diffraction,  on  the 
wavefront  curvature  induced  by  the  spherical  mirrors,  and  on  the  transverse  and  longitudinal  gain 
variations  caused  by  the  pump  mechanism.  It  would  also  be  very  desirable  to  consider  the  effects 
of  limiting  intracavity  apertures,  but  this  aspect  of  the  problem  is  more  complicated  and  has  been 
ignored  in  our  study. 

Of  course  the  mathematical  description  of  this  problem  is  considerably  more  involved  than  in 
the  case  of  a  plane-wave  model,  and  significant  theoretical  progress  is  likely  to  require  large-scale 
numerical  calculations  when  the  parameters  of  the  system  take  on  arbitrary  values.  An  unavoidable 
drawback  of  most  numerical  approaches  is  that  a  global  understanding  usually  requires  the  analysis 
of  a  wide  region  of  the  parameter  space.  For  this  reason  it  is  much  more  desirable  to  investigate 
transverse  effects  by  analytical  means,  as  much  as  feasible.  Of  course  this  cannot  be  done  for 
arbitrary  configurations  of  the  laser  system.  Fortunately  it  is  possible  to  identify  a  non-trivial 
setting  that  allows  significant  theoretical  progress  with  a  minimum  of  numerical  labor.  Our 
approach  is  based  on  an  appropriate  extension  of  the  well  known  uniform  field  limit^,  a  model 
that  has  played  a  major  role  in  defining  our  current  understanding  of  the  plane-wave  theories  of 
active  and  passive  systems.  In  spite  of  its  idealized  nature,  the  uniform  field  limit  of  the 
plane- wave  model  has  led  to  predictions  that  are  in  remarkably  close  qualitative  agreement  with  the 
results  of  the  exact  linear  stability  analysis^  and  with  the  numerical  solutions  of  the 
Maxwell-Bloch  equations^.  It  is  hoped  that  the  same  conclusion  will  remain  valid  in  the  presence 


of  transverse  effects.  Much  more  work,  however,  will  be  needed  to  confirm  this  conjecture. 

An  important  consequence  of  this  limit,  as  we  show  in  this  paper,  is  that  the  steady  state 
solutions  are  of  the  single-mode  type,  a  conclusion  which  also  implies  that  the  empty  cavity  modes 
are  also  exact  modes  of  the  filled  cavity.  This  statement  is  true  in  the  plane- wave  approximation, 
and  remains  valid  in  the  presence  of  transverse  effects.  Thus,  at  once,  the  mathematical  description 
of  the  problem  is  simplified  to  a  large  extent,  and  earlier  studies  based  on  the  single  transverse 
mode  approximation  17e,24  can  ^  understood  in  more  rigorous  terms. 

With  the  assumption  that  the  cavity  line  width  is  much  smaller  than  the  atomic  decay  rates,  the 
so  called  "good  cavity"  limit,  we  can  describe  the  emergence  of  unstable  states  through  a 
development  that  rests  almost  entirely  on  analytic  grounds;  in  fact,  only  the  final  evaluation  of  the 
rates  of  evolution  of  the  fluctuation  variables  requires  a  numerical  computation.  The  linear  stability 
analysis  leads  to  one  of  the  most  interesting  results:  the  instability  thresholds  for  a  fairly  realistic 
cavity  design  are  much  closer  to  the  first  laser  threshold  than  predicted  by  the  plane- wave  theory, 
and  are  in  much  better  quantitative  agreement  with  the  experimental  picture.  For  example,  under 
appropriate  conditions,  unstable  operation  emerges  already  when  the  pump  parameter  is  only  a  few 
percent  above  threshold,  while  the  plane-wave  theory  requires  pump  parameter  values  at  least  ten 
times  larger  than  needed  to  produce  laser  emission. 

We  show,  in  addition,  that  low  threshold  instabilities  can  be  found  also  in  the  full  adiabatic 
elimination  limit,  i.e.  when  the  dynamics  of  the  laser  is  described  entirely  by  the  field  equation 
(thus  "a  fortiori"  the  same  conclusion  holds  in  the  rate  equation  limit  when  only  the  polarization  is 
eliminated  adiabatically).  According  to  this  result  instabilities  are  not  a  necessary  consequence  of 
atomic  coherence,  in  contrast  with  the  predictions  of  the  plane-wave  theory. 

The  physical  setting  of  the  linear  stability  analysis  is  the  following:  with  the  laser  operating  in 
steady  state,  we  look  for  the  existence  of  unstable  modes,  as  evidenced  by  the  appearance  of  an 
exponential  growth  for  the  field  fluctuations.  This  approach  to  the  problem  is  conceptually  identical 
to  the  one  adopted  by  Risken  and  Nummedal^,c  in  their  classic  study  of  multimode  instabilities. 


although  there  are  also  significant  differences,  as  we  shall  discuss  in  the  main  body  of  this  report. 

In  Section  2  we  construct  the  empty  cavity  eigenfunctions  and  eigenfrequencies  for  a  folded 
ring  resonator.  In  Section  3  we  describe  the  extended  Maxwell-Bloch  model,  we  introduce  the 
uniform  field  limit,  and  calculate  the  steady  state  configuration  of  the  system  which  we  then  analyze 
for  different  values  of  the  parameters.  In  section  4  we  outline  the  details  of  the  linear  stability 
analysis.  In  Section  5  we  discuss  the  results  of  this  calculation,  give  evidence  for  the  existence  of 
low  threshold  instabilities,  and  show  how  these  effects  persist  even  after  adiabatic  elimination  of  all 
the  atomic  variables.  We  conclude,  in  Section  6,  with  a  brief  overview  and  some  general 
comments.  A  preliminary  account  of  this  work  was  published  in  Ref.  29. 

2.  The  empty  cavity  modes 

This  section  contains  a  brief  description  of  a  calculation  leading  to  the  modal  eigenfunctions 
and  eigenfrequencies  of  an  empty  resonator  in  the  paraxial  approximation^.  Its  purpose  is  to 
provide  a  self-contained  introduction  to  the  basic  tools  that  have  been  used  in  our  subsequent 
calculations. 

The  analysis  of  an  empty  resonator  of  length  A,  as  shown  in  Fig.  1,  begins  with  the 
free-field  wave  equation 

1  d 2 

V  E(r,z,t)  -  —  —  E(r,z,t)  =  0  (2.1) 

2  c  3t 

where  z  is  the  coordinate  along  the  axial  direction  of  propagation  of  the  field  and  r  is  the  radial 
coordinate  measured  from  the  axis  of  symmetry  of  the  system.  As  implied  by  Eq.  (2.1)  we 
consider  a  system  with  cylindrical  symmetry. 

In  the  paraxial  approximation  we  seek  elementary  solutions  of  the  form 

i  <V '  °V) 

E  =  Ep(r,z,t)  =  Ap(r,z)  e  p  p  (2.2) 

In  fact,  the  solutions  corresponding  to  a  cylindrically  symmetric  system  should  be  labelled  by  two 


indices  related  to  the  axial  and  radial  degrees  of  freedom,  respectively.  We  will  omit  the  axial 
modal  index,  for  simplicity,  except  where  necessary. 


In  the  slowly  varying  amplitude  approximation,  and  in  terms  of  the  scaled  coordinates 
t|  =  z/L,  p  =  /t?t  (2.3) 

v  p 

where  L  is  the  distance  between  the  curved  mirrors,  the  wave  equation  takes  the  form 

!^4(iL  +  Ii-)A  (2.4, 

dn  4  3p2  p  3p  p 

whose  solutions  are 

A  (p,T|)  =  -L  L  (■— — )  exp(-  -£?—)  exp{i  f-CL  -  (2p+l)  tan^/TlJ} 

P  V(T|)  P  v2(tl)  v2(7l)  U(T1)  ^ 

P  =  0, 1,  2, .  (2.5) 


The  origin  of  the  longitudinal  axis  is  selected  at  point  0  of  Fig.  1  where  Ep(p,Ti,t)  has  a  plane-wav 
structure.  The  functions  v(q)  and  u(T|)  are  defined  by 


v(rD  =vTl0  >/1+(T1/V 


u(Tl)=-(fi2+  fi2) 


(2.6a) 


(2.6b) 


and  T]0  =  v2(0)  s  v02  is  an  arbitrary  parameter  which  must  be  calculated  in  conformity  with  the 
cavity  geometry.  Note  that  with  the  selected  scaling  (2.3),  ri0  is  the  Rayleigh  range  of  the 
Gauss-Laguerre  beam  (2.5),  while  V  t]0  is  the  beam  waist. 

Tire  assignment  of  r)0  stems  from  the  following  considerations.  The  beam  parameter  q(r]) 
defined  by 

1  1  i 

- - - +  (2.7) 

qCn)  u(ri)  /(q) 

satisfies  elementary  mapping  rules  under  translations  and  reflections  from  the  curved  surfaces. 
Beginning  with  a  given  beam  parameter  qQ  at  position  0,  we  can  easily  construct  the  beam 


parameter  after  one  full  loop  (Fig.  1).  The  condition  q^  =  q0  secures  that  the  beam  parameter 
after  one  round  trip  acquires  the  same  starting  value. 


For  simplicity  we  select  a  cavity  with  a  high  degree  of  longitudinal  symmetry  as  shown  in 


Fig.  1 .  This  leads  to 


30 


2  1  /(  Pq-IHPq-1**) 


T1o  =  vo=2 


fd+P0)-l 


(2.8) 


where  p0  =  R/L  is  the  scaled  radius  of  curvature  of  the  spherical  mirrors  and  the  parameter  f  is 
defined  by  L/A. 

In  a  similar  way  we  can  construct  the  analytic  form  of  the  modal  functions  Ap(p,r|)  for 
l/2<T]<A/L-l/2.  The  beam  waist  in  this  case  is  given  by 


V2  = 


where  q2  is  the  beam  parameter  at  T)  =  A/2L  ( point  2  of  Fig.  1).  The  explicit  form  of  the  modal 
functions  in  this  range  is 

An(p,ri)  =  Bp[q2(n-A/2L)]  exp[-  i  (2p+l)  tan4(-  ^ L)]  e1  <p  (2.9) 


where,  for  convenience,  we  have  defined 

y^ 


2  .  /  2p2  '  '  ~2 


BB(q2(n))  - LpC -f—] )  exp(-  +  i  — £-) 

v20l)  VjfTi)  vjq)  u2(r.) 


The  functions  v2(q)  and  u2(q)  are  defined  by  Eqs.  (2.6)  with  rj0  replaced  by 

v2  -  v  W) = ■ -2-  h;  (f •  1  -  p/ + i  ( + 1  )2] 

Pq'o 


(2.10) 


(2.11) 


The  phase  factor  exp(itp)  measures  the  shift  introduced  by  the  free  propagation  along  the  segment 
1’— >2  in  Fig.  1.  This  phase  factor  can  be  calculated  by  matching  Eq.  (2.9)  at  position  1’  (i.e. 
rj= 1/2)  with  the  field  given  by  Eq.  (2.5)  after  a  reflection  from  the  first  spherical  mirror.  The 


reflection  for  near  normal  incidence  is  simulated  by  the  operator 
=  exp(-  i 

Pn 


(2.12) 


It  follows  at  once  that 


<p  =  -  (2p+l)  [tan’^— ^— )  +  tan'Vi^i-)] 
2vo  2v2 


(2.13) 


In  addition,  after  a  full  loop,  the  empty  cavity  solution  A  at  position  4  is  related  to  the  solution  at 


position  0  by  the  equation 
Ap(4)=e2i<pAp(0) 


(2.14) 


Now  we  calculate  the  empty  cavity  resonances.  The  boundary  condition  for  the  elementary 


radial  solution  or  order  p  is 

Ep(p,0,t)  =  Ep(p,A/L,x) 


i  6»nA/c 

Ap(0)  =  e  p  Ap(4) 


(2.15) 


(2.16) 


when  expressed  in  terms  of  the  modal  amplitudes  A  .  From  Eqs.  (2.16),  (2.14)  and  (2.13)  we 


A=  2nn  +  2(2p+l)  [tan  J(— ^— )  +  tan4(-^-i-)] 


(2.17) 


with  n  =  0,  ±1,  ± 2 ,....,  and  p  =,  1,  2 . The  eigenfrequencies  are  labelled  by  two  indices;  an 

axial  index  n  and  a  transverse  index  p.  The  longitudinal  intermode  spacing  is  given  by  the  usual 
formula  2nnc/A  while  the  frequency  separation  between  adjacent  transverse  modes  is 


(co  -co  )  =  5co  =4  —  [tan4(-i— J+tan'V— — -)] 

X  n,p+l  n.p'  1  .  L  v  v 


-iJ/f- 1, 


A  2v^  2m\ 


(2.18) 


An  important  property  of  the  modal  functions  Ap(p,ri)  is  that  for  every  value  of  r|  they  form 
an  orthonormal  set  in  the  sense  that 


w.y»7, 


oo 

J  dp  p  A*(p,r|)  Ap.(p,r|)  =  5p  p, 


(2.19) 


With  the  additional  reasonable  assumption  that  {Ap(p/n)}  form  a  complete  set  with  respect  to  the 
radial  coordinate,  these  functions  provide  a  useful  basis  for  the  expansion  of  any  cavity  field  at 
every  point  along  the  tj  axis.  This  property  will  play  an  essential  role  in  the  following 
development. 


3.  The  generalized  Maxwell-BJoch  problem:  description  of  the  model  and  the 
steady  state  solution 


3.1  Derivation  of  the  equations 

We  extend  the  traditional  description  of  the  laser  with  the  inclusion  of  several  features  that  are 
omitted,  by  definition,  from  the  plane-wave  theory: 


i)  the  diffraction  caused  by  the  finite  transverse  cross  section  of  the  field  and  by  its  radial 
variations  of  amplitude  and  phase; 

ii)  the  wavefront  curvature  imposed  by  the  presence  of  reflecting  spherical  surfaces; 

iii)  the  longitudinal  and  radial  variations  of  the  equilibrium  population  difference. 


We  do  so  in  the  context  of  the  paraxial  approximation  and  under  the  assumption  of  cylindrical 
symmetry  for  the  field  and  for  the  atomic  variables.  The  latter  restriction,  in  particular,  is  a 
compromise  dictated  by  the  need  for  including  transverse  degrees  of  freedom  with  the  lowest  level 
of  analytic  and  numerical  complexity.  We  recognize  that  this  is  not  a  very  realistic  setting  for  a 
typical  laser  but,  nevertheless,  we  expect  that  it  will  provide  useful  insights  into  the  role  of 


One  can  model  the  effects  of  diffraction  by  adding  a  term  of  the  form  V±F  to  the  usual  wave 
equation  for  the  slowly  varying  amplitude  F  where 


<7  52  1  3 

V  =  —  +  — ■— 

-L  -\2  r  dr 


is  the  transverse  Laplacian  for  an  axially  symmetric  geometry.  The  curved  reflectors  affect  the 
boundary  conditions  indirectly  by  way  of  the  dependence  of  the  radial  eigenfrequencies  on  the 
radius  of  curvature  of  the  mirrors  (see  Eq.  (2.17))  as  we  shall  discuss  more  explicitly  below. 
Finally,  the  longitudinal  and  transverse  variations  of  the  pump  are  reflected  in  the  equilibrium  value 
of  the  population  difference  which  is  no  longer  uniform  along  the  transverse  direction  as  in  the 
plane-wave  model. 


With  these  modifications  the  Maxwell-Bloch  equations  take  the  form 

3F  1  3F  i  ,32F  1  3Fn  ¥  „ 

3q  v  3x  4  3p2  p  3p 

^  =  -(FD  +  (l+i8AC)P) 

3x 

3D  -  r  1  *  *  1 

—  =  -Y{-t(F  p  +  fp  )  +  d-x(p,ti)} 

3t  1 


(3.2a) 


(3.2b) 


(3.2c) 


We  select  as  the  cavity  reference  frequency  co0  the  p=0  empty  cavity  resonance,  con  0,  so  that  the 
slowly  varying  field  amplitude  F  is  related  to  the  Maxwell  field  by 


E(r,z,t) 


Vll  *  ( cc  \  1  (*0 z '  “o  0  \ 

=  -  1  L  ■=•  (.F(r,z,t)  e  +  c.c.) 

Oil  z 


where  p  is  the  modulus  of  the  atomic  transition  dipole  moment  and  y1  and  Y||  are  the  polarization 
and  population  decay  rates,  respectively,  and  kg  =  (t>0/c.  The  spatial  coordinates  Tj  and  p  are  defined 
by  Eqs.  (2.3)  with  Xp  replaced  by  X0(the  reference  wavelength);  x,  v  and  8AC  are  defined  by 
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x  =  y  t,  v  =  — 

1  ly, 


SAC  = 


“a'^o 


where  ooA  is  the  atomic  transition  frequency;  the  parameter  y  denotes  the  ratio  y^/  y±  and  the  function 
X(p/n)  simulates  the  longitudinal  and  transverse  variations  of  the  pump.  There  is  a  good  deal  of 
flexibility  with  the  choice  of  x  depending  on  the  details  of  the  excitation  system.  For  definiteness 


we  select 


,  1  X(P)  M< 

x(p.n)  =  \  ZL 

l  0  otherwise 


(3.5a) 


where  LA  denotes  the  length  of  the  medium,  and  consider,  as  two  possible  models,  the  radial 


functions 


X(p)  =  exp(-  — -) 

2Pp 


X(p)  =2exp(-  -2— )-l 


(3.5b) 


(3.5c) 


In  the  first  case  the  active  medium  is  transparent  at  the  outer  boundaries  of  the  pumped  volume  (this 
model  is  plausible,  for  example,  in  the  case  of  a  four  level  system  such  as  Nd:YAG),  while  in  the 
second  a  lowering  the  pump  strength  creates  an  absorbing  region,  such  as  one  may  find  in  a  Ruby 
laser,  for  example. 

A  way  to  state  the  boundary  conditions  is  to  require  that  the  Maxwell  field  at  position  r\=0 
matches  the  field  that  has  reached  the  position  ri=A/L  after  propagating  one  full  loop  through  the 


cavity.  Formally  this  implies 

E(p,0,x)  =  E(p,A/L,x) 
or,  in  terms  of  the  slowly  varying  amplitude 


(3.6a) 


IK 

g 

» 

.T. 


F(p,0,T)  =  F(p,A/L,T)  e  u 


(3.6b) 


Equation  (3.6b)  must  be  understood  as  the  assignment  that  sets  the  slowly  varying  field  amplitude 
at  the  boundary  T|=©  and  at  time  x  in  terms  of  the  amplitude  that  has  just  completed  one  loop 
through  the  cavity  according  to  Eq.  (3.2),  and  in  compliance  with  the  appropriate  reflections  at  the 
curved  surfaces.  This  type  of  boundary  condition  is  especially  appropriate  in  connection  with  the 
numerical  solution  of  Eqs.  (3.2).  In  the  following  we  carry  out  our  analysis  of  the  steady  state  and 
of  the  linear  stability  in  terms  of  a  modal  expansion  of  the  field  so  that  a  different  version  of  the 
boundary  conditions  is  required. 

As  our  starting  point  we  let 


F(p,T|,x)  =  e  ' 80  T  F(p/n,t)  =  e ' 6n  T  A  (p/n)  f  fox) 


P(p,-n,t)  =  e'l5nx  P(p,Tl,T) 


(3.7a) 


(3.7b) 


where  8Q  is  the  offset  between  the  unknown  operating  laser  frequency  and  the  reference  frequency 
co0  in  units  of  y±.  We  substitute  Eqs.  (3.7)  into  the  Maxwell-Bloch  equations  (3.2)  and  integrate 
over  the  radial  variable  with  the  help  of  the  orthogonality  integral  (2.19).  The  result  is  the 
following  set  of  coupled  equations  for  the  modal  amplitudes  fpfox)  and  for  the  atomic  variables^ 


9fp  i  3f 


— 7  =  i  fp  -  ccLjdp  p  A*(p,q)  P(p,n,T) 


(3.8a) 


—  =  -  (FD  +  (l+iA)P) 
dx 

—  =-y{- ~-(F  p  +  FP  )  +  D  -  x(p.Tl)} 

dx  z 

where  A  =  5AC-  8Q.  The  modal  amplitudes  obey  the  boundary  conditions 


fp(-l/2,x)  =  R  e  e 


-  A-L 
iSp  i5n^i"T 


(3.8b) 


(3.8c) 


(3.8d) 


where  5p  =  (0)n  p-to0) A/c  and  R  denotes  the  reflectivity  of  the  mirrors;  a  brief  outline  of  the 
derivation  of  Eq.  (3.8d)  is  given  in  Appendix  A.  Note  that  the  sum  in  Eq.  (3.7a)  extends  to  all 
nonnegative  integers  p.  In  practice,  however,  transverse  modes  with  p  greater  than  1  or  2  suffer 
significant  diffraction  losses  typically  because  of  intracavity  apertures.  We  do  not  have  provisions 
in  our  treatment  for  describing  this  type  of  losses;  me  must  keep  in  mind,  however,  that  only  very 
few  transverse  modes  are  expected  to  play  a  significant  dynamical  role  in  the  evolution  of  a  real 
laser. 


In  steady  state  the  atomic  variables  are  given  by 


—  -  1  -i  A 

P  = - yF  , - 

st  A  si  _  2 

1  +  a2+if  I 


(3.9a) 


D  =  Y 

st  A- 


1  +  A2 


1  +  A2  +  IF 


r*  ,2 


(3.9b) 


and  the  equation  for  the  modal  amplitudes  becomes 

l  =  i42Uf  +aL(l-iA 
>  v  p 

1  0 

where 

Fs/P’7l)  =  ZAp(P’T1)fp(Tl) 


_L  =  iilfl.  f  +  aL(l-iA)  [dp  p  A*(p,rj)  x(p/n) 
dq  v  j  v 


St 


1+A2+IF  I2 

st 


(3.10) 


(3.11) 


Even  the  steady  state  problem,  in  general,  offers  formidable  analytic  difficulties.  Important 
simplifications  follow  in  the  limit 

aL. 


and 


aLA-»  0,  T  — >  0;  2C  =  -^r-  =  arbitrary 


5p-  2jim  1=0(1) 


(3.12a) 


(3.12b) 


for  all  p*0  and  m=0,  ±1,  ±2,...  The  symbol  T  =  1-R  denotes  the  mirror  transmittivity.  Equation 
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(3.12a)  is  well  known  from  earlier  studies  of  optical  instabilities  in  the  plane  wave  regime^. 
Equation  (3.12b)  implies  that  the  frequency  separation  between  adjacent  radial  modes  should  be  of 
the  same  order  of  magnitude  as  the  frequency  spacing  between  longitudinal  resonances,  and 
simultaneously  excludes  degeneracies  or  quasi-degeneracies  between  transverse  and  longitudinal 
modes.  The  two  conditions  (3.12)  taken  together  define  the  extended  uniform  field  limit  for  the 
chosen  laser  system.  Of  course,  it  is  understood  that  the  label  "uniform  field  limit"  should  not  be 
taken  literally.  What  is  actually  independent  of  rj,  in  this  limit,  is  the  set  of  modal  amplitudes  { f  } . 
The  field  profile  inside  the  cavity  is  assigned  by  the  modal  functions  { Ap)  which,  of  course, 
depend  on  the  longitudinal  coordinate  rj. 

The  importance  of  the  uniform  field  limit  in  this  context  rests  on  its  ability  to  provide  a  strong 
analytic  contribution  to  the  solution  of  the  laser  equations  of  motion  and  a  deeper  view  into  the  role 
that  the  many  parameters  play  in  setting  the  stationary  and  dynamical  properties  of  the  system. 

As  we  shall  prove  below,  the  stationary  state  of  the  laser  field  in  the  uniform  field  limit  is  of 
the  single  longitudinal  and  transverse  mode  type;  as  a  consequence,  it  follows  that  the  empty  cavity 
eigenfunctions,  in  this  limit,  are  also  modes  of  the  resonator  in  the  presence  of  the  active  medium. 

We  now  return  our  attention  to  the  solution  of  the  stationary  modal  equations  (3.10)  under 

the  conditions  (3.12).  The  spatial  integral  of  Eq.  (3.10)  yields 

+1/2  +1/2  „  _ 

f  (1/2)  -  f  (-1/2)  =  i  ML  fdn  f  +  ocL(l-iA)  fdTi  fdp  p  A*  X - =-^ - r  (3.13) 

-1/2  -1/2  0  1  +  A2  +  I FJ 

Assuming  that  5£2  is  of  the  order  of  T,  the  boundary  conditions  (3.8)  in  steady  state  can  also 
be  cast  into  the  approximate  forms 

(p=0);  f0(l/2)-f0(-l/2)=T(l-i^YL^i)  fo(i/2)  (3.14a) 


-i5„ 


A‘L 

i  sn  y  — 

1  c 


)y>  m 


(P*0);  y  1/2) '  fp<"  1/2)  =  ( 1  -R  e 


e 


(3.14b) 


J  l*|  (U  !*|  |  M  LtlJ  i ,♦ ,* *>,»  ».»  M  l.'.U* 'k*  l/.fi'-U*  •<*  la*  fc*  *,» 


If  we  combine  Eqs.  (3.14)  with  Eq.  (3.12)  we  obtain 

(P=0);  f0(l/2) - f0(-l/2)  =T(l-iMY^±!:)  f0(i/2)  = 


+  1/1  +1/1  oo  p 

i^2  fdri  f  +  aL(l-iA)  Jsq  Jdp  p  A*  %  ~ "V1 T  (3.15a) 

v  J  0  J  J  0  1+A2+IF/ 


-1/2  o 


j  5  i  5£1  — 

(p*0);  fp(l/2)  -  fp(- 1/2)  =  (l-R  e  p  e  c  )  fp(l/2)  = 


+  1/1  +1/1  ~ 
i  JQ.  J  drj  fp  +  aL(l-iA)  J  dr)  J  dp  p  A*  % 


-1/2  0 


1+A2+IF  I2 

St 


(3.15b) 


Equations  (3.15)  lead  to  the  following  conclusions: 


(i)  The  difference  f0(l/2)  -  f0(-l/2)  is  of  the  order  of  T,  so  that  f0  is  uniform  along  the  axial 
direction; 

(ii)  the  amplitudes  fp  with  p*0  is  of  the  order  of  T. 

Hence,  it  follows  that  in  the  uniform  field  limit  a  possible  steady  state  solution  is  characterized  by 
the  modal  amplitudes  f0  =  0(1),  fp*0  =  0(T)  and  the  cavity  field  configuration  is  of  the  TEM^ 


The  immediate  consequence  of  these  considerations  is  the  state  equation 


SO  A  -  L  f  f  -  'A  (p,q)l2 

7  =  2C(l-iA)  —  drj  J  dp  p  x - 2 - - 

A  JLa  o  l+A2+IA0(p,q)f0(tl)l2 


(3.16) 


for  the  modal  amplitude  fg.  where  we  have  taken  Eq.  (3.5a)  into  account.  This  can  be  split  into  real 
and  imaginary  parts  with  the  result 


ms 


oo 

=  2C^-  Jdn  Jdppx 


IA0(p.fi)l2 

i+a2+  iAfl(pji)f.(n)i2 


(3.17a) 


8£2  =  k  A  or  80  =  ■ 


(3.17b) 


where  K  =  cT/Ayx  denotes  the  cavity  linewidth  in  units  of  yx.  Equation  (3.17b)  is  an  alternative 
form  of  the  usual  mode-pulling  formula  which  remains  unchanged  in  the  uniform  field  limit  even  in 
the  presence  of  transverse  effects. 

The  stationary  state  (3.17)  is  not  unique  because,  depending  on  the  gain  parameter  2C,  other 
steady  state  configurations  are  also  possible.  In  fact  with  a  simple  extension  of  the  above 
argument,  as  shown  in  Appendix  B,  one  can  derive  the  following  more  general  result  for  a  given 
fixed  value  of  p 


=  2Cr"  J  dTi  J  dp  p  x 


1+A2+IA  f_l2 
p  p 


(3.18a) 


80  = 


ic  -  a>  5- 
=  -K-  8  +_L  .  p 

p  1+k  2rt  1+k 


(3.18b) 


where  aj  is  the  longitudinal  mode  spacing  in  units  of  yx  (aj=  2nc/A).  The  steady  state 
configuration  corresponding  to  Eqs.  (3.18)  is  also  of  the  single-mode  type;  here  the  modal 
amplitudes  f^  are  vanishingly  small  in  the  uniform  field  limit  while  f-  =  0(1).  Thus  the  solution 
of  Eqs.  (3.18)  represent  a  stationary  state  corresponding  to  the  modal  configuration  A- . 


3.2  Discussion  of  the  steady  state 

Equations  (3.18)  represent  the  main  results  of  this  section;  together  with  Eqs.  (3.9)  they 


m 


kx 


provide  a  complete  specification  of  the  steady  state  of  a  laser  and  its  dependence  on  the  many 
physical  and  geometrical  parameters  of  the  model.  We  note  that  Eqs.  (3.18)  generalize  the  results 
of  Ref.  24  and  establish  their  domain  of  validity  in  moic  rigorous  terms.  The  steady  state  formula 
derived  by  Lugiato  and  Milani  holds  in  the  uniform  field  limit  for  the  p=0  mode  and  resonant 
operation.  In  addition,  it  requires  a  flat  pump  profile  and  a  large  Fresnel  number.  To  establish  the 
connection  between  these  results  in  more  explicit  terms,  we  note  that  if  the  length  LA  of  the 
medium  is  much  smaller  than  the  Rayleigh  length  Lt|q,  Eq  (3.17a)  can  be  put  in  the  approximate 


1  =  2C 
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n  * 


exp[-  ■£—  (1  +  — ?£)] 

2p2  T1 
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,  «»/  2p2 

1  + - exp(-  — !--) 

%  % 


(3.19) 


and  if  Pp/Vq0  »  1,  Eq.  (3.19)  can  be  solved  explicitly  with  the  result 

l=2cA_in[i+l^] 

4  lfo'2 

This  agrees  with  Eq.  (4)  of  Ref.  24. 


(3.20) 


The  integrals  appearing  in  Eq.  (3.18a)  cannot  be  calculated  in  closed  analytic  form  for 
arbitrary  values  of  the  parameters.  It  is  not  difficult,  however,  to  calculate  the  threshold  gain 
parameter  for  the  p=0  mode  using  the  pump  profile  (3.5).  In  fact  after  setting  f0=0,  Eq.  (3.17) 
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(3.21) 
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Equations  (3.21)  with  =  1  and  1^=0  yield  the  threshold  gain  related  to  the  pump  profile  (3.5b); 
with  tj  =  2  and  1^=1  we  obtain  the  corresponding  result  for  the  pump  distribution  (3.5c). 

An  important  design  parameter  is  the  ratio  \{r  between  the  pump  waist  pp  and  the  minimum 
beam  size  Vrj0.  It  is  clear  by  inspection  that  the  threshold  gain  for  laser  action  decreases  as  \\i 
increases,  as  it  should  be.  The  lowest  threshold  value  corresponds  to  \jr  — >  <»  and  is  given  by  the 
formula 

TO,kt,=  (l^!)  (3.22) 

for  both  pump  models  (3.5). 

The  most  obvious  difference  between  the  two  pump  profiles  given  by  Eqs.  (3.5)  is  the 
behavior  of  the  threshold  gain  upon  varying  the  parameter  \\l.  In  the  presence  of  an  absorbing 
region  (p  >Pp)  the  threshold  gain  diverges  for  a  finite  critical  value  yc  such  that  the  denominator  in 
Eq.  (3.21 )  vanishes.  If  instead  the  active  medium  is  transparent  at  the  edges  of  the  pumped  region, 
the  threshold  gain  increases  monotonically  and  eventually  diverges  as  y  approaches  zero.  These 
considerations  are  summarized  in  Fig.  2  for  typical  values  of  r|0  and  LA/L.  [Note  that  in  order  to 
avoid  excessive  repetitions  of  the  parameters  in  the  figure  captions,  most  relevant  values  are 
tabulated  in  Table  1;  additional  comments  are  enclosed  in  each  figure  caption,  as  needed] 

The  behavior  of  the  modal  intensity  as  a  function  of  the  gain  parameter  allows  some 
interesting  observations  but  shows  no  major  surprises.  For  example,  we  looked  carefully  for  the 
possible  presence  of  bistability  near  the  threshold  for  laser  action,  as  observed  in  some  earlier 
experiments^,  but  found  no  evidence  of  this  behavior.  Some  typical  steady  state  modal  intensities 
are  shown  in  Figs.  3  and  4  for  a  resonant  and  an  off-resonant  case.  It  is  interesting  to  note  that,  for 
increasing  gain  parameter,  the  modal  intensity  of  the  radial  sidebands  can  exceed  that  of  the  p=0 
mode.  This  effect  is  related  to  the  different  transverse  profile  of  the  various  radial  eigenfunctions 


and  therefore  is  strictly  a  transverse  effect.  When  the  pumped  region  is  sufficiently  broader  than 
the  waist  of  the  TEM^  mode,  higher  order  modes  can  take  better  advantage  of  the  available  gain 
because  of  their  greater  modal  extent  (see  Fig.  3).  The  effect  is  enhanced  in  the  presence  of  a 
detuning  (Fig.  4)  because  of  the  increased  coupling  between  the  radial  sidemodes  and  the  atomic 
line.  Figures  3  and  4  correspond  to  a  pump  profile  of  the  type  given  by  Eq.  (3.5b)  with  a 
parameter  y  =  10  (i.e.  with  a  large  pump  waist  relative  to  the  beam  size).  In  this  case  one  expects 
that  the  pump  profile  (3.5c)  should  yield  rather  similar  results.  Indeed  the  differences  are  small  and 
only  quantitative  in  character. 

We  also  expect  that  the  shape  of  the  pump  profile  should  play  a  more  significant  role  when 
the  ratio  2p  /Vq0  is  closer  to  unity.  A  comparison  between  the  steady  state  curves  obtained  from 
Eq.  (3.5b)  with  those  resulting  from  Eq.  (3.5c)  is  shown  in  Fig.  5.  Here  again  the  result  is  not 
surprising  if  we  consider  that,  in  the  latter  case,  the  presence  of  an  absorbing  ring  around  the  edges 
of  the  medium  is  bound  to  affect  the  higher  order  radial  modes  more  than  it  affects  the  TEM(X) 
configuration. 

It  is  also  interesting  to  note  that  the  crossing  of  the  modal  intensity  curves  is  a  transverse 
effect  which  is  less  pronounced  when  the  intensity  profile  of  the  beam  becomes  more  homogeneous 
in  the  radial  direction,  i.e.  when  pp/Vri0»l.  An  example  is  shown  in  Fig.  6  for  a  pump  profile  of 
the  type  (3.5b). 

The  dependence  of  the  modal  intensity  on  the  detuning  parameter  8AC  is  shown  in  Fig.  7  for 
a  set  of  parameters  where  more  than  one  transverse  mode  can  be  above  threshold  for  the  same  value 
of  5ac.  In  general,  this  result  does  not  imply  that  more  than  one  mode  will  be  above  threshold  at 
the  same  time  because  the  presence  of  a  given  steady  state  alters  the  unsaturated  gain  distribution. 
Based  on  past  experience  with  the  plane  wave  model^g  we  expect  that  modes  p=0  and  p=l,  under 
the  conditions  of  the  figure,  will  compete  and  produce  either  discontinuous  jumps  (mode-hopping) 
or  undamped  pulsations  through  a  mechanism  that  was  identified  as  a  phase  instability  in  the 
plane-wave  model.  Naturally,  this  is  a  conjecture  at  this  level;  the  issue  can  only  be  verified  in  the 


context  of  the  linear  stability  analysis,  as  we  shall  discuss  in  the  next  section. 

We  conclude  this  section  by  noting  the  qualitatively  similar  behavior  of  the  output  intensity  of 
a  CO2  laser  during  a  detuning  scan,  as  shown  in  Fig.  8 

4.  Linear  stability  analysis 

The  starting  point  of  our  analysis  is  the  set  of  equations  (3.8)  for  the  modal  amplitudes  and 
for  the  atomic  variables,  together  with  the  boundary  conditions  for  the  ring  resonator.  With  the 
laser  in  a  steady  state  we  ask  whether  fluctuations  corresponding  to  the  modal  indices  (n,p)  can 
grow  exponentially  out  of  noise.  This  is  the  usual  picture  behind  the  onset  of  unstable  behavior. 
The  real  part  of  the  rate  constant  for  this  process  yields  the  rate  of  growth  of  the  unstable  modes, 
while  the  imaginary  part  measures  the  oscillation  frequency  resulting  from  the  beat  of  the  unstable 
mode  and  the  steady  state  field. 

As  our  first  step  we  map  the  coordinates  rj  and  X  into  the  new  pair  of  independent  variables 
Tp  and  x'  according  to  the  transformation 
TV  =  r| 

A-L 

T' =  x +7^  —  (q+1/2)  (4.1) 

with  the  result  that  the  new  boundary  conditions  take  on  an  isochronous  form.  Next,  we  linearize 
the  equations  of  motion  (3.8)  and  obtain 

05f  i  d5f  ^  f  *  - 

— -  + - -  =  i  —  f  -  aL  dp  p  A  (p,rj)  oP  (4.2a  ) 

dry  v  dx’  v  p  {  p 

—  =  -{5FD  ,+  F  ,5D  +  (l+iA)5P")  (4.2b) 

dx'  st  st 

~~=_v{‘  j  (tf\ + K + K + f*5*5*) + (4-2c) 


where 


5F(p,rV,x*)  =  2^  A  (p,rV)  8f  (rj',x') 


In  addition,  of  course,  one  must  also  include  linearized  equations  for  8fp*  and  8P*  which  are 
complex  conjugate  of  Eqs.  (4.2a)  and  (4.2b),  respectively.  The  boundary  conditions  are 

.  A-L 

-  i  5  ,5n^j.  ~T 

8fp(-l/2,x')  =  R  e  p  e  8fp(l/2,x')  (4.4) 

as  one  can  verify  at  once  using  Eq.  (3.8d)  in  the  new  reference  frame. 

Following  the  standard  procedure  we  introduce  the  ansatz 


SfpCn’.x1)  I 
SfJW.f) 
8P(p,rV,x') 
8P*(p,tV,x') 
18D  (p,q,,x,)J 


5<Pp(TT)  I 
ScpJW) 

eX*  8p  (p,Tj') 
8p*(p,Ti') 

8d  (p>Tj')  j 


in  the  linearized  equations  (4.2);  X  denotes  the  linearized  rate  constant  in  units  of  y±.  The  equations 
for  the  atomic  fluctuation  variables  are  algebraic  in  nature  and  can  be  solved  at  once  with  the  result 

8P(p,iy)  =  Tjtp/n’X)  a  (p/n’)  8<p  (rj  )  + 


T2(p,tT,X)  ^  A*  (p,TV)  8cp*  Cn’) 


(4.6a) 


5p*(p,n')  =  T^(p,r|',X)  ^  Ap.(p,Ti')  8cpp.(ri’)  + 

p’ 


T^pdV.X")  £  A*  (p,Ti’)  8(p*  (TV) 


(4.6b) 


where 


C3CI  "  C2C4 


CtCl  ”  C2C2 


(4.7a) 


c  jCp.Tl')  =  A.  +  1  +  iA  +  —  IFJ  (4.8a) 

1  X  +  y 

c2(p,n’)  =  Y^-Fl  (4.8b) 

z  X  +  Y 

Cjfp.n')  =  -(Ds,+  j— Fs1P*,)  (4.8c) 

X  +  y 

(4.8d> 

z  X  +  y 

The  symbols  c^*,  c-£  etc.  indicate  the  complex  conjugate  parameters  of  c^  and  c ^  etc.;  we 
note,  however,  that  in  performing  the  complex  conjugate  operation  we  must  handle  X  as  if  it  were  a 
real  variable,  i.e.  c  *  is  a  function  of  X  and  not  X*.  The  reason  for  this  peculiar  rule  can  be  traced 
back  easily  to  the  equations  of  motion  for  5p  and  5p*  and  to  the  ansatz  (4.5). 

At  this  point  the  linearized  field  equations  become 

oo 

TT6<pp  +  Yi  ^X5<pp  =  i  ^59p-aLJdpp  AplpJV)  {t^p.tV.X")  x 
Ofi  Q 

X  Ap(p’n)  5cpp(T1')  +  T2(p  ^  Ap(P’rt')  5<P*(n')  }  (4.9a) 
p'  p' 

oc 

~  5<Pp  +  Yx  7-  X  =  -  i  ^  St?*  -  cijclp  p  Ap(p,q’){T2‘(p,rT,X)x 

dTl  Q 

X  ap.(p,tv)  scpp.or)  +  t;(p,7v,x)^  a;.(P,tv)  5<p*  on }  (4.9b) 

p'  p' 


The  appropriate  boundary  conditions  take  the  form 


Before  taking  the  final  steps  with  the  calculation  of  the  linearized  eigenvalues  it  is  convenient 


to  introduce  one  more  change  of  dependent  variables  whose  main  purpose  is  to  transform  the 
boundary  conditions  to  a  standard  periodicity  form.  This  can  be  accomplished  with  the  definitions 

(il  +  t/2)  (InR  -  1 5) 


5cppW)  =  Vp(T’i')  e 


-  OV+l/2)  (InR  +  i  5  ) 


5tPp(TV)  =  Vp(TV)e 

because,  as  we  can  easily  see  with  the  help  of  Eqs.  (4.9c)  and  (4.1 1),  we  have 
vp(-  1/2)  =  vp(l/2) 

Vp(-  1/2)  =  v*(l/2) 

The  transformed  field  fluctuation  equations  take  the  form 


(4.11a) 

(4.11b) 

(4.12a) 

(4.12b) 


dv 


Ay±  drf 


oo 

(k  +  i  —  5p-i6Q)  vp+  X!  Vp  =  -  k  2C  J  dp  p  A*(pjV)  { 


Ay 


Tj(p,TT,X)  ^  Ap.(p,q')  Vp.(T| ')  e 


i  (n'+l/2)  (8p-  8p) 


T2(p,tV,X)  £  Ap.(p.Ti')  vp,(ri')  e 


i  (t|'+1/2)  (8p.+  8p)-^ 


(4.13) 


with  a  similar  equation  for  Vp*  following  from  the  complex  conjugation  of  (4.13). 

The  required  linearized  eigenvalues  cannot  be  calculated  in  closed  form  for  arbitrary  values  of 
the  parameters  in  spite  of  the  linear  nature  of  Eq.  (4. 1 3)  and  its  complex  conjugate.  However,  in 
the  good  cavity  limit,  K«l,  it  is  possible  to  develop  a  perturbative  procedure  whose  end  result  is 
the  analytic  derivation  of  the  instability  condition.  Numerical  calculations  are  needed  only  to  carry 


out  the  quadrature  of  the  final  integral  expression.  Our  strategy  can  be  summarized  as  follows:  by 
talcing  advantage  of  the  periodic  nature  of  the  boundary  conditions  (4.12),  we  expand  Vp  and  Vp*  in 
appropriate  Fourier  series  over  the  interval  -  l/2<r\<l/2  according  to  the  usual  relations 


vl' 

n 


*nLlV 
e  o 


•  >  knLll' 


(4.14a) 


(4.14b) 


where  kn  =  2rtn/L  (n  =  0,  ±1,±2,...)  and  an  and  o*n  are  complex  coefficients;  next  we 
construct  an  algebraic  system  of  equations  for  the  fluctuation  amplitudes  on  p  and  a*n  p  with  the 


help  of  the  orthonormality  relation 


1*1’ e' 


,  1  W 


(4.15) 


Finally,  we  only  need  to  solve  the  resulting  linear  algebraic  problem  to  first  order  in  K. 
Application  of  Eqs.  (4. 14)  and  (4.15)  leads  to 

i  a  o  +  (ic  +  i  — —  5  -  i  3Q)  a  +  %  o  =  k  /  ,  (  fl>n  ,p  (fc)  o  .  .  + 

n  n,p  v  A  p  n,p  n,p  >  n,p  v  '  n  ,p 


'Fn,p’(£)  <t.  .) 

nvp  v  '  n  ,p  ' 


(4.16a) 


a  a*  +  (k- i  — ^-8  +  i  5J3)  a*  +Xa*  =KT],((I)n'p  (A.)  o*.  . 

n  n,p  v  A  P  n,p  n,p  x  n*P  n  ,p 


vFn,p  \X)  a .  .) 

n,p  v  '  n  ,p  / 


(4.16b) 


where 


a  =  n  a  ,= 

n  1 


is  the  frequency  separation  between  the  n-th  longitudinal  resonance  and  the  reference  frequency 
measured  in  units  of  y±,  and  the  coefficients  and  are  defined  by 


r  kj*  x,  v  *v.  *_r 


<*<S)-2ci  ■W.,W*“’^- -V  j  dp  p  A'(p.ri')  x 


<pp'w=-2Cr-|dlVe" 


Ap.(p,Tl')  TjCp.ti'.X) 


(k„,  +  k„)LiV  i  (tV+1/2)  (5  .  +  8  ) 


(4.17a) 
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The  solution  of  Eqs.  (4.16)  is  equivalent  to  the  diagonalization  of  the  problem 


(4.17b) 


L  v  =  A.  v  (4.18) 

where  the  matrix  L,  whose  definition  is  obvious  from  Eqs.  (4.16),  has  the  structure 

L  =  Lq+kL1  (4.19) 

and  v  is  the  column  vector  with  components  an  p  and  CTn  *.  Now  we  assume  the  good  cavity 
condition  K«1  and  proceed  to  evaluate  the  eigenvalues  to  first  order  in  K.  We  assume  also  that  the 
longitudinal  mode  spacing  a|=27tc/A  and  the  atomic  linewidth  yx  are  of  the  same  order  of 
magnitude.  Hence,  the  good  cavity  condition 

k  =  (~~)t«i 


is  equivalent  to  the  uniform  field  condition  T«1 . 
We  consider  now  the  stationary  solution 

Fsi(p,rV)  =  A0(p,n’)  f0(rV) 


for  which  8f2  is  of  the  order  of  K. 


(4.20) 


(a)  Consider  first  the  case  p*0.  We  select  two  values  of  the  indices  n  and  p  (say,  n  and  pi 


t 


ft 

1 

II 


m 
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and  assume  that  there  be  no  accidental  degeneracy  of  the  type 


a  +  5  =± (a  +  — 8  ) 

B  Ay±  p  AYX 


(4.21) 


for  any  n  and  p  other  than  the  chosen  n  and  p.  Under  this  condition  the  eigenstates  of  L  have  the 


CL  _=  oL0)_  +  O (k)  =  o®  +  K  d{]_ 

n,p  n.p  n.p  n,p 


%  =  °w 

c*  -  0(k) 

n.p 


(4.22a) 


(4.22b) 

(4.22c) 


Equations  (4.22b,c)  apply  to  n,p  *  n,p.  Consider  now  Eq  (4.16a)  with  n.p  =  n.p  and  let 


X  =  X(0)  +  k  X(I)  =  -  i  (a_+  — —  8_)  +  k 

"  Ay^  P 

After  substitution  of  Eqs.  (4.22)  and  (4.23)  into  Eq.  (4.16a)  we  obtain 
-  i  (a  +-L-8ja(0)  =^°>a<0> 

n  Ay  P  n.p  n.p 


X(,)  =  -l  -V  i  SO.  +  0s.p  (^(0)) 

K  "rp 

where 


(4.23) 


(4.24a) 


(4.24b) 
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La  0 


Equation  (4.24a)  gives  the  non-degenerate  lowest  order  eigenvalues,  while  (4.24b)  gives  their  first 
order  corrections. 


(b)  Consider  now  the  case  p=0.  We  select  a  fixed  value  of  n=n  and  note  that  the  eigenvalue 


X<°>  =  -  i  a 

n 

is  doubly  degenerate,  as  we  can  easily  confirm  by  considering  Eq.  (4.16a)  for  o-0  and  Eq. 
(4. 16b)  for  a  -  0*.  It  follows  that  the  corresponding  eigenvalues  are  of  the  form 


r  =  o(0)  +  0(tc) 

n.0  n.O 

(4.25a) 
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-  “n.O  "n.O 
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(4.25b) 

r  =  O(ic) 
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for  n,p  *  n.O 

(4.25c) 

*  =0(K) 

n.P 

for  n.p  *  -n,0 

(4.25d) 

If  we  now  set  (n,p)  =  (n,0)  in  Eq.  (4.16a),  and  (n,p)  =  (-n,0)  in  Eq.  (4.16b),  and  equate  the  first 
order  terms  in  tc  we  obtain 


(i o(0)  =  <D">0a(°))  o<°)  +  4'-~n'0(xc°>)  o*  (°) 

n.O  n.O  n,0  -n.O  -n.O 
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(4.26a) 
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These  equations  lead  to  me  first  order  corrections  of  the  eigenvalues  in  the  degenerate  case.  In 
particular,  in  the  limit  La«L,T|0  and  in  resonance  (5AC=  0)  we  find  the  results  already  obtained  by 
Lugiato  and  Milani^. 

If  the  stationary  state  is  of  the  type 


Fst(p,TV)  =  Aq(p,tV)  fq(TV)  (4.27) 

correspondinr  single  radial  mode  q,  the  above  treatment  needs  only  a  few  simple 
modificatior  Equations  (4.16)  still  hold  with  replaced  by  Sfiq.The  frequency  offset  5Qq 
given  by  Eq.  t  18b)  contains  a  large  contribution  and  a  correction  of  order  k. 

The  unperturbed  eigenvalues  [Eq.  (4.24a)]  contain  only  the  large  contributions  and  are  given  by 


and  of  course  the  symbol  X(0)  that  appears  in  Eqs  (4.24b)  and  (4.26)  must  be  interpreted  according 
to  Eq.  (4.28). 

Now  we  consider  the  linearized  stability  problem  in  the  full  adiabatic  elimination  regime  of 
the  atomic  variables  (i.e.  when  y±,  y(j  »  c/A).  Here  again  the  general  procedure  needs  only  a  few 
minor  modifications  in  order  to  fit  the  limiting  situation.  In  fact,  with  reference  to  Eqs  (4.2b)  and 
(4.2c)  the  adiabatic  elimination  limit  corresponds  to  setting  both  time  derivatives  equal  to  zero.  In 
this  case  the  resulting  expressions  given  by  Eqs.  (4.8)  remain  valid  provided  we  set  X=0.  Thus  the 
linearized  eigenvalues  in  this  limit  are  still  given  by  Eqs.  (4.24b)  and  (4.26)  with  formally 
equal  to  zero  and  T  ^  and  T2  given  by  the  simple  formulas 


T/p/n)  =  -  X(P) 
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T2(P,ti)  =  X(P) 


q-iA)F;t 

(l  +  a2+  i  F  i2)2 
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In  this  case  X(1^  depends  on  p  only  by  way  of  the  modal  function  IA-12  while,  away  from  the 
adiabatic  regime,  it  depends  on  the  radial  index  also  through  the  unperturbed  frequency 


x(0) 
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P  n 
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Another  interesting  limit  is  one  in  which  the  length  LA  of  the  atomic  sample  is  much  smaller  than 
the  Rayleigh  length  Lr|0.  In  this  case  Eq.  (4.24c)  reduces  to 
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<D!,PZa(0))  =  -  2C  —  fdp  p  exp(-  He!)  L2(  )  TTp.n'.tf  01)  (4.30) 
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and  in  the  expression  for  T  j  given  by  Eqs.  (4.7a)  and  (4.8),  one  must  replace  IFst  1  with  the 
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expression 


,F/  =  ±exp(.i^)^(i^)f 


2P  'i  ,  2(  2P 


(4.31) 


It  is  important  to  remark  that  the  gain  parameter  C  in  Eq.  (4.24c)  is  related  to  the  amplitude  f- 
of  the  stationary  field,  which  appears  in  Tj  through  Fst,  by  the  state  equation  (3.18).  In  addition, 
we  note  that,  for  correctness,  the  expression  for  A  which  is  needed  in  calculating  Tj  should  be 
evaluated  to  order  zero  in  k.  Finally,  because  is  a  pure  imaginary  number,  the  instability 


condition  is  given  by  the  equation 
ReX(1>>0 


(4.32) 


5.  Stability  analysis:  Numerical  results 

The  main  results  of  the  previous  section  can  be  summarized  as  follows: 

i)  In  the  extended  uniform  field  limit  (3.12a)  and  (3.12b)  with  a,  and  yL  of  the  same  order  of 
magnitude,  a  single-mode  steady  state  field  characterized  by  modal  indices  n  and  p  becomes 
unstable  if 

Re  A/1)  >  0  (5.1) 

n.P 

for  at  least  a  pair  of  values  of  the  indices  n  and  p. 

ii)  If  p  ?t  p,  Eq.  (5.1)  takes  the  explicit  form 

-1  +ReO"’P  (X<°>)>  1  (5.2) 

where 


<&np  (X(0))  =  -  2C 

n.P  v 


oo 
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X(0)  =  -  i  (a  + - 5  ) 


(5.4) 


AY, 


and  the  function  Tj  is  defined  by  Eqs.  (4.7)  and  (4.8). 

iii)  If  p  =  p,  X(1)  is  the  solution  of  the  quadratic  equation  that  results  from  the  linear 
homogeneous  system  (4.25). 

iv)  The  instability  threshold  for  the  mode  with  indices  n  and  p,  also  known  as  the  second  laser 
threshold,  is  the  value  of  2C  that  satisfies  the  identity 


-1  +  Re  d>n-p  (X(0))=  1 
n.p 


(5.5) 


To  avoid  confusion  between  this  threshold  value  and  the  threshold  for  ordinary  laser  action,  we 
denote  the  former  by  2C<2)  and  the  latter  by  2C^. 

The  parameter  space  is  much  wider  than  in  the  case  of  the  usual  plane-wave  theory. 
Numerous  geometrical  and  physical  parameters  play  an  important  role  in  providing  the  conditions 
for  unstable  behavior  so  that  a  complete  scan  of  the  relevant  space  is  a  monumental  task.  Here,  in 
an  attempt  to  aid  our  analysis  and  to  provide  some  guidelines  for  future  searches,  we  have 
identified  the  parameters  that  can  be  viewed  as  independent  of  one  another,  at  least  from  a 
theoretical  standpoint.  For  the  general  detuned  case,  i.e.  when  ^ab  *0*  these  are  listed  in  Table  2. 
In  the  resonant  case,  and  with  the  additional  restriction  that  the  steady  state  field  is  of  the  TEM^ 
type,  the  independent  parameters  are  2C,  T)0,  \|/=Pp/Vrj0,  8o)1(  yand  LA/L;  note  that  5c0j  =  6jCtj/2rc. 

The  numerical  results  collected  in  this  section  represent  a  typical  cross  section  of  data  that  we 
have  calculated  in  a  fairly  large  size  (but  definitely  not  exhaustive)  scan  of  the  parameter  space. 
Undoubtedly,  additional  useful  facts  will  emerge  from  future  studies.  For  convenience  we  have 
divided  our  survey  into  two  subsections.  The  first  relates  to  a  resonant  laser  (SAC=  0);  the  secord 
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The  instability  threshold  is  rather  sensitive  to  the  cavity  geometry.  In  fact,  the  second  laser 
threshold  is  a  monotonic  function  of  the  transverse  mode  spacing,  as  shown  in  Fig.  9  where  we 
display  this  dependence  for  three  values  of  y.  On  the  basis  of  this  result  one  can  infer  that  long 
resonators  or  cavities  whose  curved  mirrors  have  large  radiii  of  curvature  should  be  more  sensitive 
to  instabilities  even  in  the  vicinity  of  the  threshold  for  laser  action.  In  addition,  active  media  with  a 
slow  rate  of  population  decay,  relative  to  the  decay  of  the  polarization,  are  more  unstable  than  those 
with  comparable  rates  of  decay  for  the  population  and  the  polarization.  This  behavior  is  quite 
different  from  what  is  predicted  by  the  plane-wave  theory  where  the  axial  mode  spacing  has  very 
little  influence  on  the  value  of  the  instability  threshold^. 

As  expected,  if  the  laser  operates  in  a  TEMqq  configuration  the  first  radial  sideband  (p=l)  is 
more  likely  to  become  unstable  than  the  higher  order  radial  sidemodes  (p=2,  3,...);  this  is 
illustrated  in  Fig.  10.  The  reduced  influence  of  the  higher  order  transverse  modes  on  the  unstable 
dynamics  is  linked  to  their  progressively  larger  cross  section  which  causes  the  cavity  field  to  see  a 
smaller  unsaturated  gain  because  the  equilibrium  population  difference  %  tapers  radially  out  in  a 
pre-assigned  way.  Furthermore,  in  practical  devices,  intracavity  apertures  provide  an  even  more 
stringent  selection  mechanism  by  increasing  the  losses  experienced  by  the  higher  order  modes  to 
the  point  that  they  can  no  longer  meet  the  ordinary  threshold  condition. 

One  of  the  key  parameters  that  controls  the  appearance  of  low  threshold  instabilities  is  the 
ratio  between  the  transverse  size  of  the  active  volume  and  the  beam  cross  section.  Figure  1 1  shows 
the  large  drop  of  the  instability  threshold  that  accompanies  an  increase  in  the  parameter  lor  both 
modes  p=l  and  p=2.  Figure  1  la  corresponds  to  the  pump  model  (3.5b),  while  Fig.  lib 
corresponds  to  (3.5c).  This  behavior  is  consistent  with  the  results  already  shown  in  the  previous 
figures  so  that  again,  for  example,  mode  p=2  tends  to  be  more  stable  than  mode  p=l. 

Furthermore,  the  presence  of  an  absorbing  region  at  the  radial  edges  of  the  active  medium  (Fig. 

1  lb)  causes  a  lowering  of  the  effective  ratio  between  the  transverse  size  of  the  gain  volume  and  the 
beam  waist;  hence,  the  stability  of  the  system  is  enhanced  in  the  case  shown  in  Fig.  1  lb  relative  to 
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that  of  Fig.  1  la. 


B)  Detuned  laser 

The  effect  of  detuning  is  rather  interesting  also,  although,  as  indicated  in  Table  1,  a  complete 
survey  of  the  various  possibilities  will  require  a  much  more  extensive  study  than  we  have  carried 
out  thus  far.  Figure  12  contains  a  composite  map  of  the  instability  threshold  values  plotted  as 
functions  of  the  detuning  parameter  5AC  for  four  different  off-resonant  modes,  under  the 
assumption  that  the  steady  state  of  the  laser  is  characterized  by  the  modal  indices  (n,  0).  The 
transverse  modes  (n,  1)  and  (n,  2)  exhibit  the  same  relative  stability  properties  as  displayed  in 
previous  figures,  with  mode  (n,  2)  always  being  more  stable  than  mode  (n,  1).  This  figure  shows 
also  that  transverse  modes  belonging  to  the  nearest  longitudinal  mode  (n-1)  can  play  an  important 
role  in  the  dynamics  of  this  system.  Their  resonances  are  shown  along  the  horizontal  frequency 
axis  of  the  graph  by  vertical  bars  which  indicate  that  modes  (n-1,  1)  and  (n-1,  2)  are  clustered 
around  the  resonant  mode  and  can  easily  become  unstable  as  we  vary  the  detuning  parameter. 

Mode  (n-1,  2),  in  particular,  happens  to  fall  in  the  vicinity  of  the  resonant  mode  and  has  a  low 
instability  threshold  for  a  fairly  wide  range  of  the  detuning  parameter.  We  recall,  however,  that  in 
real  laser  systems  the  losses  of  higher  order  transverse  modes  are  higher  than  those  of  lower  order 
modes,  so  that,  in  practice,  the  resonances  characterized  by  an  index  p  equal  to  or  greater  than  2  are 
likely  to  have  no  significant  dynamical  influence. 

In  an  earlier  paper  dealing  with  the  role  of  the  Gaussian  transverse  profile  on  the  dynamics  of 
a  ring  laser,  Lugiato  and  Milani^  proposed  the  surprising  result  that  the  Risken-Nummedal^'0 
instability  disappears  if  every  field  mode  has  a  TEMqq  structure.  Their  calculation  envisioned  a 
ring  laser  operating  in  resonance  and  explored  the  possible  existence  of  unstable  axial  modes  with  a 
different  value  of  the  longitudinal  index.  After  assuming  a  flat  pump  profile  (this  gives  the  best 
chance  for  the  development  of  an  instability,  as  shown  by  our  Fig.  1 1)  they  concluded  that  no 
longitudinal  sidebands  would  ever  be  unstable. 
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If  we  consider  that  the  theoretical  threshold  values  for  a  Risken-Nummedal  plane-wave 
instability  are  already  unreasonably  high,  this  result  is  very  surprising  because  it  suggests  that 
transverse  effects  have  a  stabilizing  action  on  laser  dynamics.  We  confirm  the  results  of  Ref.  24  in 
Fig.  13.  In  fact,  on  resonance  (8AC=0),  the  sidemode  (n+1, 0)  is  indeed  stable  for  all  values  of  the 
gain  (the  same  conclusion  holds  for  all  other  TEM^  modes).  It  becomes  unstable,  instead,  for 
sufficiently  large  values  of  the  detuning  parameter,  as  anticipated  in  Ref.  32,  following  a  pattern 
that  is  typical  of  the  so-called  phase  instability  of  the  plane-wave  theory  ^8.  In  the  context  of  our 
formulation  we  cannot  introduce  a  clear-cut  discrimination  between  instabilities  of  the  amplitude 
and  phase  type  because  of  the  nature  of  the  linearized  spectrum  of  eigenvalues^.  However,  the 
behavior  displayed  in  Fig.  13  is  very  reminiscent  of  the  latter  type  of  instability. 

The  combined  message  of  our  results  is  that  the  transverse  profile  of  the  field  and  of  the 
atomic  variables  does  play  an  important  role  in  the  emergence  of  low  threshold  instabilities,  but  the 
radial  modes  with  the  same  longitudinal  label  as  the  operating  steady  state  of  the  laser  are  the  ones 
that  are  most  closely  responsible  for  introducing  output  pulsations. 

6.  Conclusion  and  comments 

With  a  suitable  extension  of  the  uniform  field  limit  [Eq.  (3.12)]  we  have  succeeded  in 
carrying  out  a  detailed  analytical  study  of  the  steady  state  and  linear  stability  properties  of  a  ring 
laser  containing  a  homogeneously  broadened  active  medium.  In  this  limited  context  we  have  been 
able  to  incorporate  the  effects  of  curved  reflecting  surfaces  and  the  possible  lack  of  transverse 
uniformity  of  the  pump.  The  main  claim  of  this  paper  is  that  these  transverse  effects  can  account  ,at 
least  in  good  measure,  for  the  emergence  of  low  threshold  unstable  behavior. 

We  have  identified  the  most  critical  parameters  that  affect  the  stability  of  this  system. 
According  to  our  results,  instabilities  are  enhanced  by  a  flat  pump  profile  relative  to  the  transverse 
shape  of  the  field.  To  be  more  precise  we  found  that  if  the  ratio  2ppAJr|0  becomes  sufficiently 


large,  the  threshold  for  instability  decreases  dramatically  and  can  be  only  a  few  percent  higher  than 
the  ordinary  laser  threshold.  This  is  probably  one  of  the  most  specific  signs  of  the  transverse 
instability.  In  fact,  longitudinal  instabilities  are  favored  in  the  limiting  configuration  Pp /Vrj0 — >  0 
with  very  large  values  of  the  gain  parameter  C^e’^,  while  the  transverse  instabilities  acquire  a 
dominant  role  in  the  opposite  limit  (Pp/Vr|0— »  00  with  a  small  value  of  C).  In  this  sense,  we  can 
distinguish  between  longitudinal  and  transverse  in  stabilities.  The  former  type  is,  of  course,  a  high 
threshold  phenomenon. 

As  already  mentioned,  each  transverse  mode  with  p*0  is  associated  with  only  one 
eigenvalue,  in  contrast  with  the  situation  that  prevails  in  the  plane-wave  theory  where  two  linearized 
eigenvalues  exist  for  each  value  of  the  longitudinal  index  n.  This  precludes  a  straighforward 
distinction  between  amplitude  and  phase  instabilities.  The  obvious  cause  is  that  the  longitudinal 
modes  are  distributed  symmetrically  on  both  sides  of  the  resonant  mode,  while  the  transverse 
resonances  (p=l,2,...)  lie  only  on  one  side  of  it. 

Another  significant  difference  between  the  present  theory  and  the  its  plane-wave  counterpart 
is  the  basic  role  played  by  the  modal  functions  in  assigning  a  different  spatial  structure  to  the 
possible  field  configurations  [compare  the  Ap(p/n)  modes  with  the  functions  exp(iknLr|)].  One  of 
the  immediate  consequences  of  this  fact  is  that  all  the  field  fluctuation  variables  are  coupled  to  one 
another  to  form  an  infinite  system  of  linear  equations  [Eqs.  (4.9)].  In  contrast,  the  corresponding 
fluctuation  equations  ior  the  plane-wave  theory^?  break-up  into  an  infinite  set  of  (2x2)  blocks. 

An  important  and  fundamental  difference  between  our  model  and  earlier  plane-wave  analyses 
is  that  our  calculations  predict  the  existence  of  unstable  behaviors  in  the  rate  equation  regime 
(5P/9t=0)  and  even  in  the  limit  when  both  atomic  variables  are  eliminated  adiabatically 
(dP/5t=3D/3t=0).  This  eresult  stands  in  striking  contrast  with  the  well  known  behavior  of  the 
plane-wave  Maxwell-Bloch  equations  which  are  always  stable  in  the  rate  equation  and  in  the  full 
adiabatic  elimination  regime. 

It  is  interesting  to  compare  the  setting  and  the  results  of  this  paper  with  a  recent  analysis  by 


Lugiato  and  Lefever34  dealing  with  the  spontaneous  formation  of  space  instabilities.  An  important 
difference  is  that  the  low  threshold  temporal  pulsations  described  in  this  paper  require  a  radial 
intermode  spacing  5co1  of  the  order  of  the  longitudinal  mode  spacing  [i.e.  50),=  0(1)].  The 
emergence  of  spatial  instabilities  requires,  instead,  that  the  transverse  modes  have  overlapping 
profiles  in  frequency  and  this  implies  that  8cOj=  0(k). 

In  a  sense  this  paper  is  complementary  to  the  one  cited  in  Ref.  34;  of  course,  the  medium  is 
active  in  our  case  and  not  driven  by  an  injected  signal,  in  contrast  to  the  situation  described  in  Ref. 
34;  furthermore  the  cavity  geometry  is  of  the  ring  type  with  two  spherical  mirrors,  while  the 
resonator  of  Ref.  34  is  a  Fabry-Perot.  Both  studies,  however,  deal  with  the  competition  between 
longitudinal  and  transverse  modes.  In  our  case  the  instability  is  accompanied  by  temporal 
oscillations  with  a  frequency  of  the  order  of  the  radial  intermode  spacing;  in  the  case  of  Ref.  34,  the 
coexisting  transverse  and  longitudinal  modes  are  synchronous  with  one  another  and  produce 
stationary  spatial  patterns. 
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Appendix  A.  Derivation  of  the  boundary  conditions  (3.8d) 

With  reference  to  Fig.  1,  we  represent  the  boundary  conditions  with  the  formal  statement 

E(r,  -172,0  =  51  E(r,  172,  t  -  (A.  1) 

where  E  is  the  cavity  field,  the  coordinates  z=-L/2  and  z=L/2  correspond  to  positions  3'  and  1, 
respectively,  p  i_»3'  is  a  propagation  operator  between  the  indicated  points  and  51  is  a  reflection 
operator  that  accounts  both  for  the  reduction  in  the  field  amplitude  and  possible  curvature  effects 
upon  reflection. 

Next  we  introduce  the  slowly  varying  field  amplitude  as  in  Eq.  (3.3)  and,  with  the  help  of 
Eq.  (A.l),  we  obtain 


i<D„  A/c 

,  0  I 


F(r,  -  172,  t)  =  e  0  5t  .  F(r,  172,  t - )  (A.2) 

1— »3  C 

Because  we  are  interested  in  deriving  boundary  conditions  for  the  modal  amplitudes  fp(z,t),  we  let 
F(r,z,t)  =  e" l5<0t  A  (r,z)  fp(z,t)  (A.3) 


and,  from  Eq.  (A.2),  we  obtain 


'V~'  ico  A/c  r 

a  (r,  - 172)  f  (- 172, t)  =  e  0  e  51 


A(r,L/2)f  (L/2,t-  — )  (A.4) 

p  p  c 

The  action  of  the  reflection  operator  is  to  scale  each  modal  amplitude  f  by  an  amount  vR,  where  R 
is  the  power  reflection  coefficient,  and  to  apply  a  phase  curvature  to  the  modal  functions  A  .  If  we 
assume,  for  simplicity,  that  the  phase  mirrors  are  ideal  reflectors  and  that  the  curved  mirrors  have 
the  same  reflection  coefficient,  Eq.  (A.4)  can  be  written  in  the  form 


V  1(0  A/C  c  V  A-L 

Zi  Ap(3’}  fp(- UU)  = e  e  R  L,  fp(L/2’1  -  — } 

p  p 


SR  P^.Apd)  (A. 5) 

where  Ap(l)  and  Ap(3')  denote  the  modal  functions  at  the  indicated  points.  On  the  other  hand, 
from  the  properties  of  the  modal  functions  (See  Section  2)  we  have 


*  Ap<I)  =  V3'>e 


icop  A/c 


where  cop  is  defined  by  Eq.  (2.27),  so  that  Eq.  (A. 5)  can  also  be  expressed  in  the  form 

A-L 

v  ico  A/c  lS£1  c  Y1  A-L 

Zj  Ap(3’)  fp(-  L/2.t)  =  e  0  e  R  2,  fp(L/2,  t  -  — )  x 

p  p 

-  i(0  A/c 

Ap(3’)  e  p  (A. 

This  leads  to  the  required  result 

f  (-L/2,t)  =  Rf  (L/2,t-— )e ‘5p  (A. 

P  Pc 


where  we  have  used  the  definition  5p  =  (otp  -  co0)  A/c. 
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Appendix  B.  Derivation  of  Eqs.  (3.18) 


The  starting  point  is  provided  again  by  Eqs.  (3.10)  and  (3.1 1)  for  the  modal  amplitudes  f 
and  by  Eq.  (3.8d)  for  their  boundary  conditions.  It  is  convenient  to  introduce  a  new  set  of 
stationary  amplitude  functions  such  that 

fpCn)  =  gptn)  exp(  i  “n) 

Their  space-dependent  equation  is 


(B.l) 


9g  _  r 

— E.  =  aL(l-iA)  dp  p  A*(p,ri)  x(P,Tl) — 7  _ 

^  t  1+A2+£a_.rJ2 


A  ,  g  . 

p  bp 


(B.2) 


p 


and  the  boundary  conditions  are 


gp(-  1/2)  =  R  gp(l/2)  exp(i  (8Qy^A/c  -  5p)) 


(B.3) 


Now  we  assume  that  the  quantity 


<Pp  =  50  *  5p,  P  =  0,1,2... 


(B.4) 


is  of  the  order  of  T  for  some  integer  value  p  =  p.  In  this  case  the  boundary  conditions  (B.3)  can  be 
cast  into  the  approximate  form 


% 


(p  =  p);  g_(  1/2)  -  g_(-  1/2)  =  T  (l  -i^)  g_(I/2) 

p  p  ip 


(B.5a) 


(P  *  P);  gp(  1/2)  -  gp(-  1/2)  =  (1  -  R  e  p)gp(l/2) 


(B.5b) 


At  the  same  time  the  spatial  integral  of  Eq.  (B.2)  is 


gpCltfl-gpC-  1/2)  =  aL  (1-  iA) 


+  1/2  <» 

Jdq  Jdpp  a‘, 


ZApgP' 


-  1/2  0 


i+a2+iX 


A  .  g  .  i 

p  bP 


(B.6) 


If  we  combine  Eqs.  (B.5)  and  (B.6)  we  obtain 


(p  =  p);  g_(l/2)  -  g_(-  1/2)  =  T  (1  -i^.)g_(l/2)  =  aL(l  -  iA) 

p  p  r  p 

+1/2  «  SAp'gp' 

Jdrj  Jdp  p  A*  x  ~  -  - - 

-1/2  o  1+A2+  Ap,  gp.  I2 

p' 

(p  *  p);  gp(  1/2)  -  gp(- 1/2)  =  (1  -  R  eItPp )  ^(1/2)  =  aL  (1  -  i  A) 

+1/2  ~  SAp'gP' 

JdTi  Jdp  p  Ap  x  - - 

-i/2  o  l+A2+l2,  Ap.gp,l2 


(B.7a) 


(B.7b) 


Equations  (B.7)  lead  to  the  following  conclusions: 

i)  The  difference 

g-(  1/2)  -  g— (- 1/2) 

p  p 

is  of  the  order  of  T,  so  that  g-  is  uniform  along  the  axial  direction; 

ii)  The  amplitudes  gp  with  p*p  are  of  the  order  of  T. 

Hence,  in  the  uniform  field  limit  a  possible  steady  state  solution  is  characterized  by  the  modal 
amplitudes  gp  =  0(1)  and  g^-  =  0(T)  and  the  steady  state  is  governed  by  the  equation  (3.18) 
which  generalizes  the  case  of  the  TEMqq  steady  state  described  by  Eq.  (3.17). 
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Table  1 


Values  of  the  parameters  used  in  the  figures 


f  =  L/A 
La/L 

R  (reflectivity) 

p0  (scaled  radius  of  curvature  of  the  mirrors) 
aj  (scaled  free  spectral  range) 

5ac  (scaled  detuning) 

Vr)0  (beam  waist) 

5(0  j  (scaled  radial  frequency  spacing) 


=  0.2 

all  figures 

=  0.2 

all  figures 

=  0.9 

all  figures 

=  5.0 
=  4.1 
=  8.0 

Figs.  2,3,4,5,7,9,10,11 

Fig.  6 

Figs.  12,13 

=  1.0 
=  2.0 
=  0.5 
variable 

Figs.  2,3,4,5,6,12,13 

Fig.  7 

Fig.  1 1 

Figs.  9,10 

=  0.0 
=  0.3 
variable 

Figs.  2,3,5,6,9,10,11 

Fig.  4 

Figs.  7,12,13 

=  2.14 
=  3.36 
=  1.99 

Figs.  2,3,4,5,7,9,10,11 

Fig.  6 

Figs.  12,13 

=  0.74 
=  0.91 
=  1.48 
=  0.37 
=  0.54 
variable 

Figs.  2,3, 4,5 

Fig.  6 

Fig.  7 

Fig.  11 

Fig.  12,13 

Figs.  9,10 

=  10.0 
=  1.7 
variable 

Figs.  3,4,6,7,9,10,11,12,13 
Fig.  5 

Fig.  2 

y  =  2pp/Vq0 


Table  2 


Symbol 


Meaning 


2C 

5ac 

V=2pp/Vri0 

SWj 


Unsaturated  gain  parameter 

Difference  between  the  atomic  transition  frequency  and  the  reference 
cavity  resonance  in  units  of  y±. 

Rayleigh  range  of  the  field  profile  between  the  curved  mirrors. 

Ratio  between  the  waist  of  die  pump  medium  and  the  beam  waist  in  the 
region  between  the  curved  mirrors. 

Separation  between  the  (n,l)  and  the  (n,0)  cavity  resonances  in  units  of 

7r 


a,  Cavity  free  spectral  range. 

R  Power  reflectivity  coefficient  of  the  curved  mirrors, 

y  Ratio  between  the  population  and  polarization  relaxation  rates. 

La/L  Ratio  between  the  length  of  the  active  medium  and  the  distance  between 

the  curved  mirrors. 


Figure  Captions 


1.  Schematic  representation  of  the  ring  cavity.  The  entire  length  of  the  folded  resonator  is  A; 
the  length  of  the  active  medium  is  LA  and  the  separation  between  the  spherical  mirrors  is  L.  Point  0 
marks  the  origin  of  the  reference  system  that  measures  the  position  along  the  axis  of  the  cavity, 
while  point  4  denotes  the  end  of  a  full  loop.  Points  1  and  1\  3  and  3'  denote  the  longitudinal 
positions  just  before  and  after  reflection,  respectively.  Point  2  is  located  at  half  a  round  rip. 

2.  Ordinary  laser  threshold  gain  plotted  as  a  function  of  the  parameter  y  =  2p(/vn0  for  (a)  a 
gain  profile  of  the  type  given  by  Eq.  (3.5b)  and  (b)  for  the  gain  profile  of  Eq.  (3.5c).  As  expected, 
an  increase  in  the  transverse  dimension  of  the  excited  volume  leads  to  a  lower  threshold  gain  in 
both  cases. 

3.  Dependence  of  the  steady  state  output  intensity  on  the  gain  parameter  2C  [Eq.  (3.18a)]  for  (a) 
mode  p=0,  (b)  mode  p=l,  (c)  mode  p=2.  The  gain  profile  of  the  medium  corresponds  to  Eq. 
(3.5b).  The  larger  output  intensity  acquired  by  mode  p=l  relative  to  p=0  for  2C  greater  than  about 
2  is  a  consequence  of  the  larger  modal  volume  of  the  higher  order  modes  which  providesa  better 
match  with  the  pumped  medium. 

4.  Same  as  Fig.  3  but  with  a  detuning  parameter  5^c=  0.3  which  is  30%  of  the  chosen  free 
spectral  range  (cx  1  =  1)  and  almost  one  half  of  the  separation  between  adjacent  radial  modes 
(5coj=0.74). 


5.  Dependence  of  the  steady  state  output  intensity  on  the  gain  parameter  2C  [Eq.  (3. 1 8a)]  for  a 
medium  characterized  by  the  gain  profile  (3.5b)  for  curves  (a)  and  (b),  and  by  the  profile  (3.5c)  for 
curve  (c).  Curves  (a)  and  (b)  correspond  to  p=0  and  p=l,  respectively,  curve  (c)  corresponds  to 
p=0.  The  chosen  value  of  y  (\\f=\.7)  simulates  a  thin  column  of  active  medium.  As  expected,  in 


this  case  the  presence  of  an  absorbing  region  at  the  radial  edges  causes  a  significant  decrease  in  the 
output  intensity  of  the  p=0  mode  (the  curve  for  the  mode  p=l  is  out  of  scale). 

6.  The  dependence  of  the  output  intensity  on  the  gain  approaches  the  qualitative  behavior  of  the 
plane-wave  theory  when  the  transverse  profile  of  the  beam  becomes  more  homogeneous  in  the 
radial  direction.  Here  we  have  increased  both  the  beam  waist  and  the  width  of  the  pumped  region 
relative,  for  example,  to  the  values  used  in  Fig.  4.  Note,  however,  that  the  parameter  \j/  is  the  same 
as  in  Fig.  4  (\|/=10).  The  gain  profile  is  given  by  Eq.  (3.5b). 

7.  Dependence  of  the  modal  intensity  on  the  detuning  parameter.  Curve  (a)  shows  the  behavior 
of  the  mode  (n,p=0);  curve  (b)  corresponds  to  the  mode  (n,p=l)  and  curve  (c)  corresponds  to 
(n+l,p=0).  Note  that  the  free  spectral  range  is  0^=2  and  the  separation  between  adjacent  radial 
modes  is  5o)j  =  1.48.  The  gain  profile  is  given  by  Eq.  (3.5b). 

8.  The  output  intensity  of  a  conventional  C02  laser  as  a  function  of  the  cavity  length.  Curves 
(a)  correspond  to  the  laser  operating  in  two  consecutive  TEM^  mods  and  curve  (b)  to  a  TEM01+]0 
mode.  Additional  studies  will  be  needed  to  determine  if  the  observed  bistability  is  the  result  of 
hysteresis  in  the  piezoelectric  ceramic  driver  of  the  cavity  mirror,  or  if  it  is  intrinsic  of  the  dynamics 
of  the  system. 

9.  The  dependence  of  the  second  laser  threshold  on  the  separation  between  adjacent  radial 
modes  for  (a)  y  =  2,  (b)  y  =  1 ,  and  (c)  y  =  0. 1 .  The  gain  required  to  produce  an  instability  of  the 
p=l  mode  in  the  presence  of  a  p=0  steady  state  is  measured  in  units  of  the  ordinary  threshold  gain. 
Note  that  in  all  cases  shown  here  the  instability  threshold  is  only  slightly  higher  than  the  first  laser 
threshold  QCa)/2C(])  is  close  to  unity)  when  the  separation  between  radial  modes  becomes 
progressively  smaller.  The  gain  profile  is  given  by  Eq.  (3.5b). 
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10.  Dependence  of  the  second  laser  threshold  on  the  radial  mode  spacing  5co  j  for  (a)  p=l ,  (b) 
p=2,  and  (c)  p=3.  As  expected,  the  lowest  radial  modes  are  the  most  sensitive  to  the  appearance  of 
unstable  behavior.  The  gain  profile  is  given  by  Eq.  (3.5b). 

11.  (a)  Dependence  of  the  second  laser  threshold  on  the  parameter  y= 2pp/Vri0  for  (a)  p=l  and 
(b)  p=2.  The  pump  profile  corresponds  to  Eq.  (3.5b). 

(b)  Same  as  Fig.  1  la  with  the  pump  profile  given  by  Eq.  (3.5c).  These  results  show  that  the 
instability  is  favored  when  the  cavity  is  designed  to  operate  with  a  large  value  of  the  parameter  y. 

1 2.  Dependence  of  the  second  laser  threshold  on  the  detuning  parameter  for  a  number  of  radial 
modes;  (a)  corresponds  to  (n,  p=l),  (b)  to  (n,  p=2),  (c)  to  (n-1,  p=l),  and  (d)  to  (n-1,  p=2).  The 
pump  profile  is  given  by  Eq.  (3.5b). 

13.  Instability  domain  (to  the  right  of  the  solid  line)  for  the  radial  mode  (n+1 ,  p=0).  As  shown 
analitycally  in  Ref.  24,  this  radial  mode  cannot  become  unstable  under  resonant  conditions.  For 
sufficiently  large  values  of  5AC  and  appropriate  values  of  the  gain  parameter,  this  radial  mode 
becomes  unstable  through  a  mechanism  that  is  reminiscent  of  a  phase  instability. 
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Stimulated  Raman  scattering  from  hydrogen  molecules 


We  have  carried  out  systematic  theoretical  studies  of  stimulated  Stokes  emissions  from 
gaseous  hydrogen  molecules  irradiated  by  a  frequency-doubled  Nd:YAG  laser.  Our  objective  is  to 
investigate  the  evolution  of  the  transverse  profiles  of  the  laser  and  Stokes  fields  along  the  sample 
cell,  the  spatial  distribution  of  the  Stokes  emission,  and  the  energy  flux  exchanges  between  the  two 
interacting  fields.  This  project  was  motivated  by  some  experimental  results  obtained  by  Dr.  Y.  S. 
Liu  at  the  General  Electric  Company.  In  his  experiments  Dr.  Liu  observed  that  under  tight  focusing 
arrangements  a  single-mode  laser  produced  a  ring  pattern  of  Stokes  emission  at  a  pressure  of  60 
psi.  For  increasing  sample  pressure,  he  also  noted  the  appearance  of  a  central  peak  in  the  Stokes 
radiation,  followed  by  a  gradual  fading  of  the  ring  pattern. 


In  our  theoretical  simulation  we  assume  that  the  sample  cell  has  a  cylindrical  geometry  and 
that  the  molecular  vibration,  approximated  by  a  damped  driven  classical  harmonic  oscillator,  has 
reached  a  quasi-stationary  state  under  the  action  of  the  slower  varying  field  envelopes.  In  this  case 
the  steady-state  fields  equations  take  the  following  form: 


3FS  ia  „  g  i5  2 

—  =  — V  Fs  +  -S-(1+ - )  IF  r  F 

dfi  4  x  S  ap  1+82 


(la) 


dr  1  i  i8  2 

—  =  —  V  F  apgd - )  IFSI2  F  (lb) 

dn  4a  x  1+S2 

In  Eqs.  (1)  a  and  p  are  constants  related  to  the  ratios  of  the  frequencies  and  indices  of  refraction  of 
the  Stokes  and  laser  fields.  The  symbol  8  is  a  detuning  parameter  the  measures  the  deviation  from 
the  resonance  transition  and  g  is  a  gain  factor,  proportional  to  the  number  density  of  the  molecules. 


We  have  solved  this  set  of  equations  numerically  for  various  values  of  the  paramaters.  A 
brief  summary  of  the  results  is  given  below: 


1.  For  a  large  input  beam  waist  only  a  central  peak  is  obtained  in  the  Stokes  emission  at 
the  exit  of  the  sample  cell.  In  this  case,  generally,  the  spatial  pattern  of  the  transmitted  laser  field 
shows  both  a  central  peak  and  a  ring-shaped  structure.  For  a  smaller  input  beam  waist  a  ring 
structure  appears  in  the  spatial  distribution  of  the  Stokes  radiation  and  the  energy  transfer  is  so 
complete  that  the  laser  field  is  totally  depleted. 

2.  If  the  Stokes  wavelength  is  in  resonance  or  shifted  to  the  red  side  of  the  molecular 
transition,  no  ring  pattern  develops  in  the  Stokes  emission.  This  feature  becomes  evident  only 
when  the  Stokes  radiation  is  blue-shifted  relative  to  the  molecular  transition. 


3.  In  order  to  simulate  the  effect  of  varying  pressure  we  have  investigated  the  dependence 
of  the  solutions  on  the  gain  parameter.  For  appropriate  selection  of  the  remaining  parameters  of 
the  system,  a  ring  pattern  appears  for  small  gain;  for  larger  gain  values  the  central  peak  begins  to 
develop  and  eventually  becomes  a  dominant  feature  of  the  spatial  distribution  of  the  Stokes 
emission.  These  results,  as  well  as  item  1,  are  consistent  with  Liu's  experimental  observations. 

4.  The  ring  pattern  in  the  Stokes  emission  may  disappear  when  the  input  laser  intensity 
becomes  sufficiently  small.  On  the  other  hand  this  pattern  is  insensitive  to  the  weak  level  of  the 
Stokes  intensity  assumed  at  the  entrance  of  the  sample  cell. 
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Miscellaneous  Theoretical  Results 

1 .  We  have  completed  a  theoretical  and  experimental  analysis  of  the  behavior  of 
inhomogeneosly  broadened  ring  lasers.  The  results  of  this  investigation  include  the  quantitative 
analysis  of  the  behavior  of  the  self-pulsing  frequency  as  a  function  of  the  most  important  control 
parameters  and  a  study  of  the  pulsation  waveforms.  In  a  paper  published  in  the  Physical  Review  A 
by  Tarroja,  Abraham,  Bandy  and  Narducci  we  have  shown  that  heterodyne  and  homodyne  spectra 
together  with  phase-spae  plots  can  be  useful  in  categorizing  the  periodic  and  chaotic  behavior  of 
such  systems.  The  theoretical  model,  developed  in  earlier  studies  in  collaboration  with  Lugiato  of 
the  Liniversity  of  Milano  (presently  at  the  Polytechnic  Institute  of  Torino)  is  successful  in 
explaining  a  large  number  of  experimental  observations  carried  out  in  our  laboratories. 

2.  Albano  has  implemented  new  techniques  for  calculating  dimensions  and  entropies  from 
digitized  data.  There  is  increasing  evidence  that  broad-band  laser  spectra  may  be  a  symptom  of 
dynamical  chaos,  while  similarly  broad  spectra  in  amplified  spontaneous  emission  may  only  be  the 
reflection  of  stochastic  noise. 

3.  Yuan  has  completed  a  study  of  the  time  evolution  equation  for  the  density  matrix  of  a  driven 
Morse  oscillator  coupled  to  a  bath  of  harmonic  oscillators.  The  results  show  that  the  expected 
hysteresis,  which  is  typical  of  classical  calculations,  does  not  appear  in  the  context  of  quantum 
mechanical  calculations.  On  the  other  hand,  bistability  is  manifested  in  the  bimodal  population 
distribution  attained  by  the  nonlinear  oscillator.  Distributions  of  this  type  have  ben  observed  by 
Bloembergen  and  others.  Time  series  of  the  oscillator  energy  show  increasingly  more  complicated 
behavior  as  the  driving  intensity  increases. 

4.  New  results  have  been  obtained  by  Yuan  in  the  study  of  the  quantum  behavior  of  a 
dynamical  system  which  is  known  to  exhibit  classical  chaos.  In  particular,  we  have  studied  the 
dissociation  dynamics  of  a  nondissipative,  laser-driven  Morse  oscillator  as  a  prototype  of  IR 
multiphoton  dissociation  of  molecular  systems.  In  the  classical  treatment  of  this  problem  we  found 
primary  resonances  with  winding  numbers  equal  to  1/n  (where  n  is  a  small  integer),  corresponding 
to  the  fundamental  and  overtone  transitions  of  the  Morse  oscillator.  Furthermore,  an  infinite 
sequence  of  secondary  resonances  are  generated  as  a  result  of  the  interaction  of  1 :1  and  1:2  primary 
resonances.  We  have  found  also  that  the  dissociation  half  line  near  the  threshold  of  global 
instability  follows  a  power  law  relation,  with  the  critical  exponents  in  close  agreement  with  those 
predicted  by  the  standard  map.  This  fact,  the  dissociation  rate  and  and  the  dissociation  fractions  are 
amenable  to  a  direct  comparison  with  experimental  and  quantum  theoretical  resul's. 

5.  We  have  used  calculations  of  the  dimension  of  strange  attractors,  correlation  functions  and 
higher  order  coherence  properties  to  analyze  and  identify  the  nature  of  the  dynamical  or  stochastic 
causes  of  intensity  fluctuations  in  lasers  and  amplified  spontaneous  emission  sources.  We  have 
obtained  experimental  data  from  long  amplifiers  and  unidirectional  ring  lasers.  A  preliminary 
account  of  these  techniques  was  included  in  M.H.F.  Tarroja's  PhD  Dissertation. 

6.  We  have  developed  and  analyzed  a  model  for  bidirectional  homogeneously  broadened  ring 
laser  such  as  should  describe  our  previous  experiments  in  CO2  lasers.  We  have  found  new  results 
for  time-dependent  periodic  and  chaotic  solutions  and  we  have  identified  the  boundaries  for  stable 
operation.  We  have  also  extended  the  traditional  models  of  the  single  mode  homogeneously 
broadened  ring  laser  to  new  operating  conditions  and  found  new  characteristic  behaviors.  We  have 
analyzed  a  model  of  a  laser  with  saturable  absorption  under  conditions  in  which  the  absorber  and 
amplifier  media  have  different  resonant  frequencies  and  when  the  cavity  is  detuned  from  the  center 
frequency  of  the  active  atoms.  The  results  of  our  calculations  predict  both  bistable  and  unstable 
solutions. 
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7.  We  have  implemented  a  modification  of  the  Grassberger-Procaccia  algorithm  for  calculating 
the  dimensions  of  a  strange  attractor  using  the  singular  value  decomposition  techniques  of 
Broomhead  and  King.  This  new  method  is  especially  useful  in  verifying  if  a  dynamic  attractor  is 
truly  fractal  over  all  scale  lengths. 


Experimental  Results 


1 .  Following  a  theoretical  suggestion  by  P.  Mandel  that  a  laser  swept  from  below  to  above 
threshold  would  display  discontinuous  jumps  and  hysteresis,  Tredicce  and  Scharpf  have  produced 
an  experimental  demonstration  of  this  effect  with  an  Ar+  laser.  Because  the  original  theoretical 
proposal  was  phrased  in  such  terms  that  a  quantitative  comparison  was  not  readily  possible, 
Tredicce  developed  a  new  theoretical  description  of  this  process  and  obtained  good  agreement 
between  the  experiments  and  the  results  of  his  calculations.  These  investigations  have  been 
published  in  Optics  Communications.  Lippi,  Hoffer  and  Abraham,  at  Bryn  Mawr,  are 
investigating  the  same  effects  in  He-Xe  and  He-Ne  ring  and  Fabry-Perot  lasers. 

2.  By  controlling  the  level  of  losses  inside  tha  cavity  of  an  Ar+  with  the  help  of  an  acousto-optic 
modulator,  Tredicce  has  studied  the  spectral  behavior  of  the  laser  output.  The  result  of  an  analysis 
carried  out  with  a  Fabry-Perot  interferometer  has  shown  that,  as  a  function  of  decreasing  losses,  at 
frrst  the  laser  operates  in  a  single  mode,  then  a  second  mode  appears,  and  finally  a  third  mode  also 
tries  to  exceed  its  threshold  condition.  At  this  point,  usually,  chaotic-like  behavior  sets  in.  Further 
investigations  are  in  progress  to  clarify  the  nature  of  the  broad-band  spectral  features. 

3.  The  Drexel  Quantum  Optics  Group  has  completed  construction  of  two  C02  lasers.  The  first 
operates  with  a  power  output  of  about  15  Watts,  and  is  used  mainly  to  study  transverse  effects  in 
laser  instabilities.  The  second  laser  is  a  low  power  (1-2  Watts)  sealed-off  model  which  operates  as 
a  reference  source  for  heterodyne  measurements. 

4.  We  have  carried  out  research  to  seek  improved  characterization  of  both  ring  and  Fabry-Perot 
laser  behaviors.  We  have  measured  power  spectra  of  output  intensity  and  digitized  the  intensity 
pulsations  using  the  spectrum  analyzer  and  Lecroy  fast  transient  digitizer  which  have  been 
purchased  under  the  subcontract  to  Bryn  Mawr.  Detailed  studies  of  the  digitized  data  have  provided 
new  ways  to  measure  chaos  using  correlation  function^  and  intensity  statistics. 

5.  Gioggia,  Abraham  and  their  consultant  Lange  (from  the  University  of  Hanover,  FRG)  have 
analyzed  the  Fabry-Perot  laser  searching  for  the  recently  predicted  regions  of  bistability  and 
instability  in  the  laser  frequency  near  the  Lamb-dip  region.  Considerable  progress  has  been  made  in 
understanding  the  data.  A  paper  containing  these  results  will  be  published  in  the  special  issue  of 
the  Optical  Society  of  America  devoted  to  laser  instabilities. 
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Special  Activities 


1.  Abraham  Co-Chaired  a  Symposium  on  Chaos  which  was  held  in  Quebec  as  part  of  the  1986 
SPIE  Meeting.  The  Symposium  ran  for  two  days  and  featured  a  program  of  invited  and  contributed 
talks  from  practically  all  the  major  contributors  in  the  area  of  laser  instabilities  [Arecchi,  Abraham, 
Boyd,  Casperson,  Chrostowski  (the  second  Co-Chairman),  Lugiato,  Harrison,  Narducci, 

Tredicce,  Weiss,  and  several  others].  Because  many  of  the  participants  were  also  involved  in  the 
joint  CLEO/IQEC  meeting  in  San  Francisco,  this  Symposium  provided  the  opportunity  to  Fine  tune 
some  of  the  presentations  that  were  scheduled  to  be  given  to  a  much  broader  audience  within  a  few 
days. 

2.  In  the  Quebec  Symposium  on  Chaos  Abraham  gave  an  overview  of  ways  to  characterize 
chaotic  phenomena  in  lasers,  and  discussed  new  ways  to  measure  the  dimensionality  of  an  attractor 
and  to  characterize  chaotic  behavior  in  dynamical  systems. 

3.  In  the  same  Symposium  Lugiato  and  Narducci  proposed  a  new  classification  scheme  for 
labelling  the  unstable  behaviors  that  arise  in  laser  and  laser-related  systems  based  on  the  concepts  of 
amplitude  and  phase  instability.  Their  work  is  limited  to  the  plane-wave  approximation,  but  it 
represents  a  significant  improvement  relative  to  earlier  approaches  based  on  the  concepts  of  single 
and  multimode  operation. 

4.  Tredicce  discussed  a  number  of  peculiar  instabilities  observed  in  COj  lasers  in  collaboration 

>  with  his  Florence  colleagues.  He  attributed  the  appearance  of  these  pulsations  to  transverse  effects 

and  proposed  a  phenomenological  description  based  on  the  concept  of  intensity-dependent  cavity 
damping  rate;  the  qualitative  agreement  of  the  resulting  theoretical  analysis  and  the  experimental 
results  is  quite  remarkable.  The  theoretical  analyses  leading  to  these  results  were  carried  out  jointly 
by  the  Bryn  Mawr  and  Drexel  groups. 

5.  Abraham,  Mandel,  and  Narducci  have  just  completed  a  full  length  review  paper  on  Laser 
Instabilities.  This  article  is  in  press  and  will  appear  as  part  of  the  series  Progress  in  Optics.  The 
manuscript,  420  typewritten  pages  long,  consists  of  10  chapters  and  more  than  800  references. 

6.  A  volume  of  lecture  notes  on  Laser  Physics  and  Instabilities,  by  Narducci  and  Abraham,  is 
in  press  by  the  World  Scientific  Publishers.  This  work  is  based  on  the  lectures  given  by  the  authors 
at  a  Summer  School  in  Laser  Physics  held  in  Changchun,  China,  in  August  1985. 

7.  Abraham,  Boyd,  and  Arimondo  organized  a  Conference  entitled  Optical  Instabilities  II  which 
was  held  in  Italy  in  July  1987.  Narducci  was  a  member  of  the  program  committee. 

8.  Narducci  was  elected  Fellow  of  the  American  Physical  Society  for  his  theoretical 
contributions  to  Quantum  Optics  and  especially  optical  instabilities. 

9.  J.R.  Tredicce  and  D.K.  Bandy  are  co-editors,  with  A.  Oraevskii  of  the  Lebedev  Institute  of 
Moskow,  of  a  special  issue  of  the  Journal  of  the  Optical  Society  of  America  devoted  to  recent 
progress  in  the  nonlinear  behavior  of  laser  systems. 

10.  J.M.  Yuan  was  the  organizer  of  the  "Chaos  in  Molecular  Systems"  Symposium  of  the  '86 
International  Laser  Science  Conference  held  in  Seattle,  Washington,  Oct.  1986. 


1 1.  J.M.  Yuan  was  one  of  the  organizers  of  the  XV  International  Colloqium  on  Group 
Theoretical  Methods  in  Physics  held  at  Drexel  University,  Oct.  1986. 
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12.  N.B,  Abraham,  L.M.  Narducci  and  J.M.  Yuan  gave  invited  talks  at  the  International 
Workshop  Optical  Instabilities  II  held  at  II  Ciocco,  Italy  in  July,  1987. 

13.  N.B.  Abraham  and  L.M.  Narducci  were  among  the  lecturers  at  a  NATO  Summer  Institute  on 
Laser  Physics. 

14.  N.B.  Abraham  was  promoted  to  the  rank  of  Full  Profesor  at  Bryn  Mawr  College. 

15.  L.M.  Narducci  was  nominated  Francis  K.  Davis  Professor  of  Physics  and  Atmospheric 
Science  at  Drexel  University 

16.  L.M.  Narducci  has  become  a  member  of  the  Center  for  Advanced  Studies  at  the  University 
of  New  Mexico,  Albuquerque. 

17.  J.M.  Yuan  was  an  invited  speaker  at  the  2nd  International  Laser  Science  Conference,  held  in 
Seattle,  Washington,  Oct.,  1986. 

18.  J.M.  Yuan  was  an  invited  speaker  at  the  Telluride  Summer  Research  Institute  on  Nonlinear 
Dynamics,  June  1987. 

19.  J.M.  Yuan  was  a  Symposium  Leader  at  the  Sanibel  Symposium  held  at  Marineland,  Fla., 
March,  1987. 
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